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ABSTRACT 



Spitzer IRAC observations of two fields in the extended UV disk (XUV-disk) 
of M83 have been recently obtained, ~3 Rhii away from the center of the galaxy 
{Rhii=Q.6 kpc). GALEX UV images have shown the two fields to host in-situ 
recent star-formation. The IRAC images are used in conjunction with GALEX 
data and new HI imaging from The HI Nearby Galaxy Survey (THINGS) to 
constrain stellar masses and ages of the UV clumps in the fields, and to relate 
the local recent star formation to the reservoir of available gas. Multi-wavelength 
photometry in the UV and mid-IR bands of 136 UV clumps (spatial resolution 
>220 pc) identified in the two target fields, together with model fitting of the 
stellar UV-MIR spectral energy distributions (SED), suggest that the clumps 
cover a range of ages between a few Myr and >1 Gyr with a median value 
around < 100 Myr, and have masses in the range 10^-3x10^ M©, with a peak 
~ 10^'^Mq. The range of observed ages, for which only a small fraction of 
the mass in stars appears to have formed in the past ~10 Myr, agrees with the 
dearth of Ha emission observed in these outer fields. At the location of our IRAC 
fields, the HI map shows localized enhancement and clumping of atomic gas. A 
comparison of the observed star formation with the gas reservoir shows that the 
UV clumps follow the Schmidt-Kennicutt scaling law of star formation, and that 
star formation is occurring in regions with gas densities at approximately (within 
a factor of a few) the critical density value derived according to the Toomre Q 
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gravitational stability criterion. The significant 8 /im excess in several of the 
clumps (16% of the total by number accounting for ~67% of the 8 /xm fiux)) 
provides evidence for the existence o f dust in these remote fields, in agreement 
with results for other galaxies (e.g. iPopescu fc Tuffj l2003l ). Furthermore, we 



observe a relatively small excess of emission at 4.5 /im in the clumps (14%±6% 
by fiux), which suggests contribution from hot small grains (~ lOOOK), as already 
observed in other galaxies. From our data, the outer regions of the M83 galaxy 
disk show evidence of a time-extended star formation history over <1 Gyr, and of 
a moderately chemically-evolved interstellar medium, in agree ment with recent 



findin gs on the metallicity of the outer HII regions of M83 (iGil de Paz et al. 
2007h . 



Subject headings: galaxies: evolution 
ISM — galaxies: star clusters 



galaxies: individual (M83) — galaxies: 



1. Introduction 



Star formation in the outer disks of galaxies has multiple implications for our under- 
standing of the formation and evolution of disks, of the laws that govern star formation, 
and of the interaction of massive stars with the rarefied interstellar medium in those fields. 
Presence of star formation ensures that those external fields are undergoing some chemical 
enrichment. Radiative and mechanical feedback from massive stars may be more efficient in 
a low-density environment, which may thus play a key role in the enrichment of the pristine 



halo. The star formation proce ss in the outer disk takes place late compared to the inner 

Furthermore, the low gas density enables tests of the star 



disk (iMuhoz-Mateos et al 



formation threshold (IMartin fc Kennicuttll200ll ) and of the relation betweeri gas density and 
star formation rate density (the Schmidt-Kennicutt Law, iKennicutt I Il998l ) at the low end 
of the range. Overall, outer disks provide insights into low-density conditions for the star 
formation that may have characterized the early disk formation. 

Deep Ha imaging had already revealed outer disk star formation beyon d two optical 
radii (-R25) in a few nearby galaxies (IFerguson et al.lll998l : iLelievre fc Royll2000l ). Broadband 
obser vations had also shown presence of sig nificant numbers of B stars in the outer disk of 
M31 JCuillandre et allboOlh and NGC 6822 Jde Blok et al.ll2003h . Unlike their counterparts 
in the inner disks, the outer disk HII regions are small, faint and isolated. Thus, the presence 
of suc h HII regions does not repre sent a challenge to the notion that Ha 'edges' exist in most 
disks dMartin fc Kennicuttl[200ll ). 
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More recently, GALEX has shown that the UV pr ofile of disks extends, with a smooth 



trend, well beyond the Ha radius or the R25 of galaxies (IThilker et al.ll2005l : iGil de Paz et al. 



20051 . Thilker et al., 2007, in press, Boissier et al., 2007, in press). XUV-disks are apparently 
common in the present epoch - though not ubiquitous -, as t hey are found at least at a 
low level in about one-quarter of the spiral galaxy population (IZaritsky fc ChristleinI 120071 . 
Thilker et al. 2007, in press). The UV clumps in the outer disk are generally associated with 
extended HI struct ures, and show eyidenc e for metal enrichment, with metallicities in the 



range Z^/S-Zq/IO flGil de 



within extended HI disks fiThilker et al 



'az et al 



2007). Indeed, extended UV disks are generally hosted 



20051 : iGil de Paz et al.ll2005l ). Thus, the rarified 



outer disk is not quiet and star formation process is occurring stealthily in the extended HI 
disk. The presence of smooth UV surface brightness profiles does not counter the paucity of 
Ha emission found in the outer disk UV emission, since the UV probes a larger range of ages 
than Ha and at low-star- formation levels the number of ionizing stars may be very small 
and stochastically absent (IBoissier et al.l 120071 ). Thus, the outer disk star formation merits 
careful investigation, to understand its nature and the conditions under which it occurs. 

The southern grand-design galaxy M83 (RA(2000): 13h37m00.9s, Decf2000): -29d5 1m57s) 



is a metal rich (> Z q) Sc galaxy, viewed almost face on ( i^ 25", ISofue et al.lll999l ). Its 
distance, 4.5 Mpc (IThim et al.l l2003l : iKarachentsevI 120051 ). makes the spatial resolution 
of Spitzer (2", ~44 pc) and GALEX (5", ~109 pc), sufficient to isolate large star for- 
mation complexes. The mean star formation rate per unit area of M83 is ro ughly 0.04 



Mf7\ yr ^ kpc ^ and the total star formation rate (SFR) is about 5 Mq yr ^ (IKennicutt 



19981 ). Therefore, M83 is relatively active among the local normal star forming galaxies. 
The Ha emission of M83 has a wel l-defined 'edge' at the galactocentric distance of around 
6.6 kpc (IMartin fc Kennicuttll200ll ). M83 also has a very extended HI disk which is about 
3 times larger than the optical one (R25). Th e possibihty of a clos e tidal encounter with its 
companion NGC 5253, around 1-2 Gyr ago (IRogstad et al.lll974l : Ivan den Bergh.l Il980l ) is 
still controversial, but if confirmed, it could have triggered the starbursts in both galaxies. 
Additio nal sources of interaction f or this galaxy include dwarf companion galaxies such as 
KK208 jKarachentseva et al.lE998[ ). 



GALEX images of M83 have revealed the presence of more than 100 UV-bright sources 
grouped in highly str uctured complexes beyond the radius where the Ha surface brightness 
decreases abruptly (IThilker et al.ll2005l ). However, only a few of these UV bright regions 
(10%-20%) have Ha counterparts. This apparent discrepancy has been sugges ted to be due 



to th e wide range of the stellar populations' age within the outer disk of M83 (IThilker et al. 



20051 ) and/or to the stochastic nature of the stellar IMF at the very low gas densities in 
these remote regions (IBoissier et al.l 120071 ). 
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New mid-infrared images, from Spitzer IRAC , of two fields in the extreme outskirts of 
M83 are presented here, and are combined with t he GALEX iraages and a 2 1 cm (HI) image 
from The HI Nearby Galaxy Survey (THINGS, IWalter et all J2005I . boosh ) to investigate 
the properties of the UV-bright extended disk; we are exploiting the leverage offered by the 
multi-wavelength observations to constrain ages and masses, and to infer the star formation 
characteristics of the UV knots in the outer regions of M83. The targeted IRAC fields 
show relatively high UV surface brightness and are associated with large scale filamentary 
HI structures. The THINGS HI image offers about 3 times better angular resolution than 
pre-existing HI images, and enables us to explore in detail the scaling laws of star formation 
at a local level. 



In section 2, we provide a brief description of our data. We describe in section 3 the 
methods used to select the sample of regions for the detailed analysis. The analysis of the 
photometry is presented in section 4 and our results are summarized in section 5. 



Observations and Data Reduction 



Spitzer images at 3.6, 4.5, 5.8 and 8.0 /im of two UV bright f ields in the outsk irts of 
M83 were obtained during 2005, July using the IRAC instrument (IFazio et al.ll2004l ). The 
two fields, called Outerl and OuterL, are selected for the combination of both UV-brightness 
and Ha deficiency (see Fig. [1]) and detailed information of these fields is listed in Table [1] 
Spitzer IRAC observed each field twice with a time interval of several days to enable removal 
of asteroids. The two fields were mapped with 4x7 and 4x6 grids, for Outerl and OuterL, 
respectively, with sub-frame steps in order to reach maximum depth in fields of size of 
7.5' X 7.5' and 7.5' x 5'. The total observation time per pixel in the center of each map is 480 
second, reaching a depth of 0.007 MJy/sr at 3.6 /im (1 a detection error). The standard 
Spitzer Infrare d Nearby Galaxies Survey (SINGS) IRAC pipeline was employed to create the 
final mosaics (IKennicutt et al.ll2003l : iRegan et al.ll2004l ). 



M83 was observed by GALEX in 2003 in both the far-UV (FUV, A ~15 29 A) and the 
near-UV (NUV, A ;^231 2 A), as part of the Nearby Galaxy Survery (NGS; iBianchi et al. 



20031 : iMartin et al.ll2005l ). The total exposure time is 1352 s in each filter and the Point 
Spread Function has a FWHM=4.6^^ The detailed description of the GALEX observations 
is provided in iThilker et al.l (120051 ). The GALEX images were aligned, registered, and 
cropped to match the Field-of- Views (FOVs) of the IRAC images of Outerl and OuterL, 
respectively, (Fig. [1]) and re-sampled to the same pixel size (0.75" /pixel for Outerl and 
1" /pixel for OuterL). The GALEX and Spitzer images of Outerl and OuterL are depicted 
in Fig. [2] and Fig. [31 respectively . 
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The 21 cm HI image of M83, from THINGS flWalter et all l2005l . l2008h . was reduced 
and calibrated according to standard pro cedures developed for the THINGS project (e.g. 
Walter et al.l l2005l : iKennicutt et al.l 120071 ). with the addition of a careful blanking in the 
data cube to improve sensitivity in the outer regions of the map. Furthermore, the map has 
been corrected for the primary beam attenuation, as the two IRAC fields lay at the edges of 
the half-power point of VLA primary beam (~30' diameter). The final beamsize is 15"xll". 
The final HI image has also been aligned, registered, and cropped to match both the Outerl 
and Outer! FOVs. 



Images from the 2MASS survey (ISkrutskie et al.ll2006l ) at J, H, and were also used 
in the present analysis to provide a sanity check on our results. The 2MASS images are at 
least 2 orders of magnitude less deep than our IRAC images, and the resulting uncertainty 
on the photometry of our sources (which are very faint in the 2MASS fields) prevents the 
effective use of the 2MASS data in our fitting routines (see below). However, the 2MASS 
data can still provide consistency checks, and they will be used as such. The 2MASS images 
corresponding to the locations of our fields where retrieved and aligned to both Outerl and 
OuterL. 



3. Source Identification and Selection 

Outerl and OuterL contain a few hundred UV sources, which can be, in addition to 
star forming regions in M83, background galaxies and foreground stars. Identification of the 
nature of the sources and of the "bona fide" IRAC counterparts to their UV emission (rather 
than chance superpositions) has required extensive investigation. We developed procedures 
to both exclude distant galaxies and foreground stars (often showing in both the UV and mid- 
IR, M83 lies at an low ecliptic latitude of ~ 18°), and flag UV sources that are contaminated 



Table 1. Spitzer IRAC observations 



Target 


Field Center Position 


Galactocent. 


IRAC Map's 


exposure 


Field 


(J2000) 


Distance 


core size 


tinie(s)'^ 


M83 Outerl 


13:35:56.02 -29:57:23.0 


15.06' (19.71 kpc=2.96 Rhii) 


7.5'x7.5' 


3000 


M83 OuterL 


13:36:56.90 -30:06:41.7 


14.77' (19.33 kpc=2.91 Rhii) 


7.5'x5.0' 


2573 



''For the central map locations, the effective exposure is 480 second 
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in the IRAC bands by chance superposition of physically unrelated mid-IR sources. We 
applied several steps, as discussed below, to achieve these goals. 

We began constructing our sample by selecting FUV- and NUV-emitting regions from 
the GALEX images in order to include all potential star forming regions. The FUV is sen- 
sitive to the young phases of SF (though not effectively instantanenous like Ha) , while 
the NUV traces the slightly evolved stell ar population (age ~ 100 Myr). Sextractor (ver- 



sion2.4.4, iBertin. E.. fc Arnouts. S.lll996l ) was employed to identify all 2a sources, yielding 
226 and 133 in Outerl and OuterL, respectively (after manually combining several adjacent 
and overlapping apertures in the crowded regions). Since the instrument background and 
sky background of GALEX and Spitzer are much different, we did not use the Sextractor 
photometry package to measure the flux of these apertures in the IRAC bands. Instead, 
the IRAF routine "phot" was used to extract the photometry of the UV-detected sources in 
all the six GALEX+ IRAC bands after subtracting the local background from each of them. 
The radius of each region was selected >5" (~the GALEX PSF, corresponding to a physical 
scale ~220 pc in diameter). The same position and aperture diameter was used in all six 
(FUV, NUV, 3.6, 4.5, 5.8 and 8.0 iim) bands and aperture corrections are applied to the re- 
sulting photometry. All our sources appear point-like at the resolution of both GALEX and 
Spitzer/IRAC; therefore we apply to each source/bandpass combination the point-source 
aperture correction appropriate for the chosen region size. For the GALEX measurements, 
we derived the aperture corrections from photometry with increasing aperture radii of point 
sources in the M83 field (outside the galaxy); as a sani t y che ck, we compared our results with 



the aperture corrections published by lMorrissey et al.l (120071 ). and found no major discrepan- 



cies with our values. For the IRAC measurements, we referred to the point-source aperture 



corrections published in the IRAC Data HandbooHj. Some sources are dim in several bands, 
and after subtracting the local background, the flux is negative. For these la upper 

limit has been set. 

Uncertainties on the measurements for each source/bandpass combination are a quadratic 
combination of three contributions: variance of the local background, photometric calibration 
uncertainties, and variations from potential mis-registration of the multiwavelength images. 
We find that in all cases the variance of the local background is the predominant source of 
uncertainty for our measurements, so we only consider this contribution to the uncertainties 
quoted in Table [2l The local background is derived from fits of the pixel distribution of 
regions of a few thousands pixels (or about l/50th of the area of each image) immediately 



5IRAC Data Handbook, v3.0, Ch. 5; 



http: / / ssc.spitzer.caltech.edu/irac/dh/ 
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surrounding each aperture, for which we derive the mode and the variancqj. Each region 
was also visually inspected to ensure that no systematic trends in background areas were 
present. In addition, tests were run on sample apertures to verify that the mode of each local 
background is not dependent on the size of the area selected for the background measure- 
ments: we varied region sizes by up a factor of 2 (both increasing and decreasing) without 
encountering for the vast majority of cases variations on excess of the measured variances. 

Contaminating foreground and background sources are eliminated through several meth- 
ods. For some nearby stars, we can identify them by their shapes in the 3.6 /im images, i.e., 
by the presence of diffraction spikes. We also exploit HST ACS images of the center parts 
of the two fields, obtained by one of us (D.T.), to further identify and remove contaminated 
sources from the sample; the ACS images, thanks to their superior angular resolution, en- 
able a quick recognition of foreground stars and background galaxies. For the most part, we 
drop from our candidates list any source with strong contamination. In those cases where 
the contaminating source is physically much smaller than our measurement aperture (e.g., 
a distant background galaxy) and close to its edge, we shift and resize the measurement 
aperture to avoid it. In a few instances, a foreground star may lay close to the center of our 
selected aperture: in these cases, we assume that the FUV and NUV fluxes receive negligible 
contamination, and we only modify our IRAC photometric measurements, by subtracting 
the contaminating contribution measured in an inset small radius aperture. In all cases, we 
apply the point-source aperture correction appropriate for the aperture size selected. 

Three especially bright MW stars lie within the FOVs of Outerl (to the right) and 
OuterL (both right and left), see Fig. [2] and Fig. [31 They have strong wings due to the PSF 
of the Spitzer IRAC instrument. We reject all apertures that fall within and close to these 
wings. After all identified contaminating sources are removed from or corrected for in our 
initial list, there are 166 and 114 regions remaining in the Outerl and OuterL, respectively. 

Finally, in order to remove dim foreground low mass stars more completely, we adopt 
two strong selection criteria: log{^^) > —0.5 and fpuv > 5o", where fpuv and Jnuvi in 
unit of ergs cm~^ A~^, are the flux of FUV and NUV separately and a is the measurement 
uncertainty. Although the criteria will cause the loss of some clusters with age larger than 
~ 1 Gyr (from the prediction of Starburst99, for an instantaneous burst of star formation 
with 1/5 solar abundance), it is very effective at eliminating low mass stars which emit most 
of their energy in the optical and near-Infrared band. This 'loss' of the oldest clusters will 



^The routine used for the fits, the IRAF task MSKY written by M. Dickinson (1993, private commu- 
nication), allows the user to interactively define the interval in the pixel distribution where the mode and 
variance are calculated. This ensures that robust results are obtained even in the absence of source masking. 
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have minimal impact on this analysis, which centers on more recent star formation events 
(younger than about 1 Gyr). 

The final size of our sample is 152 sources: 97 in Outer! and 55 in OuterL. The final 
source list with their fiuxes and 1 a uncertainties is listed in table [2l Although we try our 
best to exclude the possible foreground stars and background galaxies, there are still about 
16 sources showing unusual high excess 4.5, 5.8 and 8.0 yum fiux which cannot be easily 
explained by stellar population synthesis models (see section These sources tend to 

be among the bright ones in our sample, and we cannot attribute the observed excess to 
photometric or background subtraction uncertainties. We separate these unusual sources 
from the other normal sources in table [2l We consider the remaining 136 sources in our final 
list (although we will carry along the other 16 sources from some of the comparisons, for 
completeness.) 



Table 2. Source List 







RA(J2000) 


Dec(J2000) 




FUV'' 


NUV'' 


3.6^ 


4.5^ 


5.8^ 


8.0'' 


ID 


IRAC field 


degree 


degree 


Radius* 


Fx (lo-) 




Fx 


Fx (!<') 


Fx (l") 


Fx 



1 


Outerl 


13 


36 


22 


1 


-29 


52 


22 


1 


11 


4U4cs.4( yy 


2) 


2700 


2 


38 


5) 


45 


6 





65) 


17 


64 





65) 


7 


27 


1 


96) 


5 


2 


0.91) 


2 


Outcrl 


13 


36 


14 


3 


-29 


53 


05 


5 


11 


2814.4(106 


3) 


1198 


2 


41 


2) 


39 


5 





70) 


19 


79 





70) 


9 


85 


2 


09) 


19 


3 


0.97) 


3 


Outcrl 


13 


36 


19 


2 


-29 


53 


4t) 


8 


10 


I361 1 .Z{ \)Z 


2) 


7391 


8 


35 


8) 


15 


1 





60) 


9 


48 





60) 


11 


75 


1 


81) 


1 1 


3 


0.85) 


4 


Outerl 


13 


36 


26 


5 


-29 


54 


10 


5 


8 


id * 4. ( ( ( 


9) 


779 


7 


30 


2) 


1 1 


4 





51) 


7 


72 





51) 


2 


75 


1 


54) 


5 


1 


0.72) 


5 


Outerl 


13 


36 


17 





-29 


54 


52 


1 


8 


y&o.y^ fi 


0) 


418 


5 


27 


4) 


17 


3 





47) 


8 


27 





47) 


2 


67 


1 


4t)) 


3 





t).65) 


6 


Outerl 


13 


35 


56 


6 


-29 


54 


54 


1 


14 


7945.2(127 


5) 


3041 


1 


49 


3) 


31 


2 





84) 


13 


60 





84) 


4 


76 


2 


51) 


6 


1 


1.17) 


7 


Outerl 


13 


36 


07 


5 


-29 


55 


15 


3 


11 


3810.2(106 


3) 


2926 


3 


41 


2) 


156 


6 





70) 


71 


91 





70) 


24 


98 


2 


09) 


29 


4 


0.97) 


8 


Outerl 


13 


36 


12 


2 


-29 


57 


06 





8 


2987. 6( 71 


0) 


1245 


5 


27 


4) 


57 


2 





47) 


25 


08 





47) 


11 


65 


1 


40) 


2 


8 


0.65) 


9 


Outcrl 


13 


36 


08 


9 


-30 


01 


11 


3 


8 


1249. 1( 71 


0) 


887 


7 


27 


4) 


20 








47) 


9 


11 





47) 


6 


05 


1 


40) 


2 


7 


0.65) 


10 


Outcrl 


13 


35 


49 


4 


-30 


t)3 


32 


5 


8 


I/J5.4( i L 


0) 


978 


8 


27 


4) 


80 


1 


{) 


47) 


40 


95 





47) 


11 


70 


1 


40) 


8 


8 


t).65) 


11 


Outerl 


13 


35 


51 


1 


-30 


00 


16 


3 


8 


1095. 4( 71 


0) 


1015 


9 


27 


4) 


117 


8 





47) 


49 


38 





47) 


22 


27 


1 


40) 


7 


6 


0.65) 


12 


Outerl 


13 


35 


48 


4 


-29 


58 


48 


8 


8 


1030. 5( 70 


8) 


486 





27 


4) 


30 


9 





47) 


15 


73 





47) 


3 


88 


1 


39) 


5 


6 


0.65) 


13 


Outerl 


13 


35 


25 


5 


-29 


58 


13 





8 


878. 9( 70 


8) 


759 


4 


27 


4) 


21 


6 





47) 


7 


74 





47) 


3 


68 


1 


39) 


9 


6 


0.65) 


14 


Outerl 


13 


35 


40 


5 


-29 


58 


31 


4 


8 


iiyo.o( fi 


0) 


826 


9 


27 


4) 


41 


6 





47) 


24 


10 





47) 


8 


08 


1 


40) 


14 


3 


0.65) 


15 


Outerl 


13 


35 


46 


5 


-29 


53 


37 


8 


11 


2857.7(106 


3) 


1765 


2 


41 


2) 


58 


4 





70) 


31 


35 





70) 


11 


35 


2 


09) 


18 


2 


0.97) 


16 


Outcrl 


13 


35 


35 


8 


-29 


51 


24 


8 


6 


874. 6( 56 


9) 


489 


4 


21 


9) 


24 


4 





37) 


11 


72 





37) 


4 


58 


1 


12) 


15 


4 


0.52) 


17 


Outcrl 


13 


36 


27 




-29 


54 


30 


8 


8 


588. 9( 71 


0) 


968 


7 


27 


4) 


171 


6 





47) 


74 


66 





47) 


28 


05 


1 


40) 


1 1 


2 


0.65) 


18 


Outerl 


13 


36 


06 


1 


-29 


57 


55 


5 


6 


759. 9( 56 


7) 


803 


3 


21 


9) 


19 


1 





37) 


11 


11 





37) 


3 


54 


1 


11) 


10 





0.52) 


19 


Outerl 


13 


36 


09 


3 


-29 


56 


38 


2 


9 


876. 8( 85 


1) 


1036 


2 


32 


9) 


55 


7 





56) 


31 


35 





56) 


11 


40 


1 


68) 


11 


2 


0.78) 


20 


Outerl 


13 


35 


45 





-30 


00 


07 


9 


6 


1026. 2( 56 


7) 


641 


3 


21 


9) 


27 


5 





37) 


14 


54 





37) 


4 


87 


1 


11) 


9 


6 


0.52) 


21 


Outcrl 


13 


35 


39 




-30 


01 


34 


8 


6 


617. 0( 56 


7) 


567 





21 


9) 


14 


8 





37) 


9 


19 





37) 


3 


06 


1 


11) 


7 


3 


0.52) 


22 


Outcrl 


13 


35 


50 




-30 


02 


39 


6 


8 


326. 9 ( 71 


0) 


742 


6 


27 


4) 


12 


3 





47) 


5 


04 


t) 


47) 


2 


75 


1 


4t)) 


1 


9 


t).65) 


23 


Outcrl 


13 


35 


37 


2 


-29 


51 


15 


8 


8 


439. 5( 71 


0) 


810 


1 


27 


4) 


243 


8 


t) 


47) 


95 


82 





47) 


35 


00 


1 


40) 


1 1 


2 


0.65) 


24 


Outcrl 


13 


36 


04 


5 


-29 


53 


10 


7 


8 


1125.7( 7t) 


8) 


637 


9 


27 


4) 


2 


8 





47) 





87 





47) 


2 


83 


1 


39) 





6 


0.65) 


25 


Outerl 


13 


36 


04 


3 


-29 


54 


32 


5 


11 


3680.3(106 


3) 


1586 


4 


41 


2) 





6 





70) 





44 





70) 


2 


69 


2 


09) 





7 


0.97) 


26 


Outerl 


13 


36 


16 


5 


-29 


55 


33 


7 


8 


1309. 8( 71 


0) 


398 


3 


27 


4) 


12 


5 





47) 


3 


68 





47) 


2 


36 


1 


40) 





4 


0.65) 


27 


Outerl 


13 


36 


25 


6 


-29 


57 


01 


2 


9 


1680. 0( 85 


1) 


418 


5 


32 


9) 


82 








56) 


33 


51 





56) 


14 


49 


1 


68) 


1 


5 


0.78) 


28 


Outerl 


13 


36 


20 


3 


-29 


56 


46 


5 


8 


1428. 8( 71 


0) 


381 


4 


27 


4) 


4 


2 





47) 


3 


31 





47) 





68 


1 


40) 





6 


0.65) 


29 


Outcrl 


13 


36 


02 


3 


-29 


56 


15 


3 


8 


1980. 9( 71 


0) 


789 


8 


27 


4) 











47) 





32 





47) 





85 


1 


40) 





3 


0.65) 


30 


Outcrl 


13 


35 


59 


1 


-29 


56 


23 


5 


6 


2084. 8( 56 


7) 


594 





21 


9) 





2 





37) 





05 





37) 





71 


1 


11) 


1 


3 


0.52) 


31 


Outerl 


13 


35 


58 


7 


-29 


57 


10 


4 


6 


A A O 1 Ci^ 

^y44.d( OD 


9) 


1340 





21 


9) 





8 





37) 


1 


14 





37) 


1 


53 


1 


12) 





6 


0.52) 


32 


Outerl 


13 


35 


57 


7 


-29 


56 


31 


3 


6 


iOyy.Ot^ Do 


9) 


489 


4 


21 


9) 





3 





37) 





26 





37) 





71 


1 


12) 





5 


0.52) 


33 


Outerl 


13 


35 


56 


6 


-29 


56 


37 


5 


11 


13465.7(106 


3) 


6075 


4 


41 


2) 


8 


1 





70) 


5 


94 





70) 


3 


33 


2 


09) 


2 


9 


0.97) 


34 


Outerl 


13 


35 


55 


3 


-29 


56 


40 


3 


6 


2294. 8( 56 


7) 


850 


6 


21 


9) 


1 


2 





32) 





73 





32) 





42 





97) 





4 


0.45) 


35 


Outcrl 


13 


36 


06 


5 


-29 


59 


48 


4 


9 


2468. 0( 85 


1) 


718 


9 


32 


9) 





3 





56) 





25 





56) 





97 


1 


68) 





6 


0.78) 


36 


Outcrl 


13 


35 


56 


3 


-29 


58 


10 


7 


8 


1801. 2( 70 


8) 


745 


9 


27 


4) 


4 


3 





47) 


1 


55 





47) 





33 


1 


39) 





3 


0.65) 


37 


Outcrl 


13 


35 


54 


4 


-29 


58 


34 


3 


8 


Vi-Zl 1 L 


0) 


391 




27 


4) 


2 


2 





47) 


3 


06 





47) 


1 


00 


1 


40) 





2 


0.65) 


38 


Outcrl 


13 


35 


55 





-29 


59 


51 


2 


11 


4199.9(106 


3) 


1518 


9 


41 


2) 


8 


3 





70) 


3 


59 





70) 


1 


54 


2 


09) 





1 


0.97) 


39 


Outerl 


13 


36 


10 


5 


-30 


00 


38 


3 


8 


989. 4( 71 


0) 


739 


2 


27 


4) 


4 


7 





47) 


2 


41 





47) 


1 


98 


1 


40) 


1 





0.65) 


40 


Outerl 


13 


36 


15 


8 


-30 


01 


54 


9 


8 


1260. 0( 71 


0) 


607 


5 


27 


4) 


6 


9 





47) 


4 


21 





47) 


1 


93 


1 


40) 


7 





0.65) 


41 


Outerl 


13 


36 


04 


4 


-30 


01 


57 


Q 




1 ^f\Ck QQ 


2^ 


570 




38 


5^ 


10 


Q 


Q 


65) 




66 


Q 


65) 




93 




96) 


Q 




91) 


42 


Outerl 


13 


36 


01 





-30 


01 


43 


1 


8 


1253. 5( 71 


0) 


546 


8 


27 


4) 


15 


3 





47) 


6 


84 





47) 


1 


85 


1 


40) 





4 


0.65) 


43 


Outcrl 


13 


35 


57 


9 


-30 


01 


50 


2 


11 


5065.9(106 


3) 


2420 





41 


2) 


19 


7 





70) 


8 


21 





70) 


2 


90 


2 


09) 


1 


8 


0.97) 


44 


Outerl 


13 


36 


00 


7 


-30 


04 


04 


5 


11 


3290.6(106 


3) 


1758 


5 


41 


2) 










70) 





95 





70) 





75 


2 


09) 





3 


0.97) 


45 


Outerl 


13 


35 


36 





-29 


59 


11 


8 


8 


840. 0( 70 


8) 


421 


9 


27 


4) 


9 


1 





47) 


3 


35 





47) 





57 


1 


39) 


2 


1 


0.65) 


46 


Outerl 


13 


35 


44 


9 


-29 


58 


03 


4 


8 


1608.5( 71 


0) 


567 





27 


4) 


3 


6 





47) 


1 


11 





47) 





33 


1 


40) 





1 


0.65) 


47 


Outerl 


13 


36 


14 


6 


-29 


53 


56 


9 


8 


346. 4( 71 


0) 


749 


3 


27 


4) 


16 








47) 


8 


43 





47) 


2 


09 


1 


40) 


8 


6 


0.65) 


48 


Outerl 


13 


36 


15 


3 


-29 


52 


48 


7 


8 


378. 9( 71 


0) 


843 


8 


27 


4) 


25 


6 





47) 


10 


93 





47) 





66 


1 


40) 


1 


7 


0.65) 


49 


Outcrl 


13 


36 


010. 


-29 


52 


37 


4 


5 


287.9( 42 


4) 


418 


5 


16 


5) 


6 


4 





28) 


2 


82 





28) 





27 





84) 





5 


0.39) 


50 


Outcrl 


13 


35 


58 


6 


-29 


52 


27 


5 


5 


242. 5( 42 


4) 


325 


4 


16 


5) 


3 


1 





28) 





91 





28) 





36 





84) 





1 


0.39) 


51 


Outerl 


13 


35 


40 


6 


-29 


51 


30 


9 


5 


456. 8{ 65 


2) 


1026 


1 


24 


1) 


23 


3 





41) 


8 


37 





42) 





91 


1 


22) 








0.59) 


52 


Outerl 


13 


36 


07 


6 


-29 


58 


40 


5 


5 


257.6( 42 


4) 


371 


3 


16 


5) 


32 


1 





28) 


12 


66 





28) 


5 


39 





84) 





4 


0.39) 


53 


Outerl 


13 


36 


07 


9 


-30 


00 


15 





5 


521. 7( 49 


6) 


658 


2 


19 


2) 


4 


4 





33) 


2 


33 





33) 


1 


96 





98) 





4 


0.45) 


54 


Outerl 


13 


35 


52 


1 


-29 


59 


42 


2 


5 


519. 6( 42 


4) 


408 


4 


16 


5) 


1 


3 





28) 





28 





28) 


1 


06 





84) 





2 


0.39) 


55 


Outerl 


13 


36 


07 


1 


-29 


51 


41 


1 


8 


592. 6( 71 


0) 


358 


8 


27 


4) 





9 





47) 





89 





47) 





85 


1 


40) 





6 


0.65) 


56 


Outerl 


13 


36 


14 


6 


-29 


54 


31 


5 


15 


3788.6(141 


7) 


1066 


6 


54 


8) 


16 


2 





93) 


12 


61 





93) 


3 


40 


2 


79) 


2 


3 


1.30) 


57 


Outcrl 


13 


36 


24 


8 


-29 


58 


00 


1 


8 


708. 8( 71 


0) 


315 


4 


27 


4) 


3 


8 





47) 





58 





47) 


3 


75 


1 


40) 


3 


4 


0.65) 


58 


Outerl 


13 


35 


27 


9 


-30 


01 


36 





8 


1027.9( 71 


0) 


331 


6 


27 


4) 


70 


2 





47) 


27 


10 





47) 


11 


14 


1 


40) 


2 


1 


0.65) 


59 


Outerl 


13 


35 


42 


3 


-29 


53 


56 


5 


8 


790. 9( 71 


0) 


514 


1 


27 


4) 


11 


1 





47) 


6 


94 





47) 


1 


66 


1 


40) 


5 


1 


0.65) 


60 


Outerl 


13 


36 


20 


4 


-29 


52 


05 


9 


6 


338. 5( 56 


7) 


281 


9 


21 


9) 


4 


3 





37) 


3 


69 





37) 


4 


58 


1 


12) 





8 


0.52) 


61 


Outerl 


13 


36 


22 


1 


-29 


54 


30 


7 


10 


1530.4( 92 


1) 


805 


6 


35 


7) 


8 


2 





61) 


2 


85 





60) 


1 


02 


1 


81) 





7 


0.84) 


62 


Outerl 


13 


36 


21 


1 


-29 


56 


03 


8 


8 


1021. 1( 71 


0) 


291 


9 


27 


4) 


3 








47) 





37 





47) 


1 


70 


1 


40) 





8 


0.65) 


63 


Outcrl 


13 


36 


11 


1 


-29 


56 


10 


5 


8 


893. 4( 71 


0) 


521 


4 


27 


4) 


13 


3 





47) 


6 


85 





47) 


1 


70 


1 


40) 


3 


8 


0.65) 


64 


Outcrl 


13 


36 


01 


8 


-30 


00 


04 


5 


8 


647. 3( 71 


0) 


333 


1 


27 


4) 


6 


2 





47) 


1 


86 





47) 





85 


1 


40) 


2 


1 


0.65) 


65 


Outerl 


13 


35 


36 


6 


-30 


02 


28 


1 


8 


806. 8{ 71 


0) 


211 


1 


27 


4) 


13 


7 





47) 


6 


35 





47) 


1 


35 


1 


40) 


2 


8 


0.65) 



- 10 - 



I 




Fig. 1. — [adapted from iThilker et al.l (120051 )]: The GALEX image of M83 in shown in a 
two-color combination, blue for the F UV channel and orange for the NUV channel. The red 
contours outline the extent ed HI disk fiTilanus fc Allenlll993l ). The yellow contour indicates 
'^neutral gas = 10 Mq pc~^ ( Crosthwaite et al. 2OO2I ). The two green squares at the right and 
the bottom of picture indicate the positions and approximate size of our two Spitzer IRAC 
fields; the name convention for the two fields is also indicated on the figure. North is up, 
East is left. 
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Fig. 2.— GALEX (FUV, NUV) and Spitzer IRAC (3.6 and 8.0 //m) images of the Outer! 
field. The black circles show the 97 sources retained in the sample after removal of contami- 
nating foreground and background sources. The size of the circles indicates the apertures in 
which photometry is performed. 




Fig. 3. — As in Figure [21 now for the OuterL field. A total of 55 sources are retained for 
this field. 
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Table 2 — Continued 



RA(J2000) Dec(J2000) FUV'' NUV*" 3.6^ 4.5'' B.s'' S-O*" 

ID IHAC field degree degree Radius* (Ict) (Ict) F^ (la) F^ (Ict) Fx (Ict) F^ (Ict) 
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4. Analysis 



4.1. Stellar populations' modelling 

The age and mass of each source is constrained by fitting the UV+MIR bands with 
models of the spectral energy distribution (SEDs) of single-age stellar populations. This 
implicitly assumes that the emission in the IRAC bands is due to pho tospheric eniission 



from stars. This is a reasonable assumption for the 3.6 and 4.5 /im bands (jPahre et al.l 12004 



Calzetti et al.ll2005l ). The 5.8 and 8 fim band potentially contain a significant (or dominant) 
contribution from dust emission, due to Policyclic Aromatic Hydrocarbons single-photon 
heated by the stellar populations; to evaluate the impact of dust emission on our stellar 
population fits, we create a dust-only 8 /xm image by rescahng and subtracting the 3.6 /im 
image from the 8 fim one. The resulting dust-only image shows a number of sources with 
detected dust emission, in most cases a weak detection, in agreement wi th the expectation 



that low-metallicity regio ns have underluminous 8 fim dust emission ( lEngelbracht et al. 



20051 : ICalzetti et al.l 120071 ). However, to prevent out fitting results from being even mildly 



contaminated by dust emission contributions, we use only the two IRAC bands at the shortest 
wavelengths in our fits, together with the two GALEX bands (next section). 

The theoretical SEDs are frorn the Starburst99 stellar population synthesis models 
( jLeitherer et al.lll999l : IVdzquez et al.ll2005l ). We select instantaneous burst populations, since 
we are likely to deal with single or small groups of stellar clusters (i.e., same or similar age 
populations) due to the small physical sizes of the regions sampled (lOOpc to 300pc). A 
metalhcity value of Z=0.004 is adopted for the models, to match as closely as possible the 
oxy gen abundances of roug hly 1/10-1/5 Zq measured in a few of the HIl regions in these ar- 
eas (iGil de Paz et al.ll2007l ). The Kroupa Initial Mass Function ($( m)(xm~°, with a=1 .3 for 
O.IMq < M < 0.5Mq, and a=2.3 for 0.5Mo < M < lOOM©, i JKroupa et allboOlh . the 
default IMF of Starburst99, is employed. Since we include in the fits the mid- Infrared (MIR) 
fluxes of each stellar complex, we use the new Padova stellar models with full asymptotic 
giant branch evolution. In order to balance the accuracy and the speed of the simulations for 
different age ranges, we select different time steps for the models: 0.05 Myr for 0.1-1 Myr, 
0.1 Myr for 1-10 Myr, 0.5 Myr for 10-50 Myr, 1 Myr for 50-200 Myr, 2 Myr for 200-500 Myr, 
10 Myr for 500-1000 Myr and 50 Myr for 1-2 Gyr. We convolve the spectral response curve 
of the GALEX and Spitzer filter bandpasses with the simulated stellar population spectra 
to obtain the synthetic photometry in the six bands. Finally, the luminosities and the colors 
from observations are compared with the theoretical model to constrain ages and masses. 
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Table 2 — Continued 







RA{J2000) 


Dec(J2000) 




FUV'' 


NUV'' 


3.6'' 


4.5'' 


5.S^ 


8.0*= 


ID 


IRAC field 


degree 


degree 


Radius^ 


Fx (l"^) 




Fx {1<^) 


Fx (1<^) 


Fx (l"^) 


Fx (1<t) 


131 


OuterL 


13 37 14.8 


-30 09 


55.7 


7 


773.0(118.7) 


546.5(33.4) 


5.1(0.68) 


0.67(0.66) 


0.92(1.82) 


0.3(1.04) 


132 


OuterL 


13 36 56.5 


-30 07 31.0 


6 


1022.4(101.8) 


630.9(28.6) 


0.4(0.58) 


0.21(0.57) 


0.67(1.56) 


0.3(0.89) 


133 


OuterL 


13 36 58.0 


-30 07 


52.9 


5 


823. 0( 84.8) 


509.5(23.9) 


1.2(0.48) 


0.56(0.47) 


0.41(1.30) 


0.5(0.74) 


134 


OuterL 


13 36 57.9 


-30 08 


20.8 


5 


635. 2( 84.8) 


486.8(23.9) 


2.8(0.48) 


1.22(0.47) 


1.44(1.30) 


0.4(0.74) 


135 


OuterL 


13 36 57.1 


-30 07 


59.0 


6 


536.5(101.8) 


873.6(28.6) 


65.8(0.58) 


27.69(0.57) 


9.00(1.56) 


2.5(0.89) 


136 


OuterL 


13 36 57.1 


-30 09 


19.8 


5 


665. 3( 84.8) 


352.5(23.9) 


3.5(0.48) 


1.81(0.47) 


0.89(1.30) 


0.8(0.74) 










Sources with unusual high excess 


4.5, 5.8 and 8 


flux 








137 


Outerl 


13 36 03.4 


-29 55 


25.7 


8 


1420. 2( 70.8) 


742.6(27.4) 


57.5(0.47) 


29.98(0.47) 


15.08(1.39) 


39.5(0.65) 


138 


Outerl 


13 35 55.1 


-29 56 


57.5 


7 


3702. 0{ 63.9) 


1937.4(24.7) 


13.7(0.42) 


7.54(0.42) 


8.18(1.26) 


5.7(0.58) 


139 


Outerl 


13 36 13.4 


-29 59 


42.1 


8 


1474. 3{ 71.0) 


2028.5(27.4) 


61.2(0.47) 


61.33(0.47) 


51.61(1.40) 


35.7(0.65) 


140 


Outerl 


13 35 49.9 


-30 00 


40.7 


8 


5022. 6{ 71.0) 


5299.1(27.4) 


56.0(0.47) 


55.45(0.47) 


50.67(1.40) 


40.0(0.65) 


141 


Outerl 


13 35 28.4 


-29 58 


39.7 


8 


1907.3{ 71.0) 


1437.9(27.4) 


52.0(0.47) 


21.95(0.47) 


13.67(1.40) 


24.9(0.65) 


142 


Outerl 


13 35 39.8 


-29 58 


48.3 


8 


1405. 0{ 71.0) 


1066.6(27.4) 


75.6(0.47) 


43.50(0.47) 


19.21(1.40) 


65.1(0.65) 


143 


Outerl 


13 36 15.6 


-29 58 


00.0 


8 


493. 6{ 71.0) 


864.1(27.4) 


93.3(0.47) 


60.74(0.47) 


41.36(1.40) 


68.2(0.65) 


144 


Outerl 


13 35 48.9 


-29 56 


40.9 


8 


532. 6{ 71.0) 


617.7(27.4) 


24.3(0.47) 


20.77(0.47) 


15.79(1-40) 


9.4(0.65) 


145 


Outerl 


13 35 34.1 


-29 52 


25.5 


8 


703. 6{ 71.0) 


621.0(27.4) 


61.5(0.47) 


34.68(0.47) 


18.38(1-40) 


20.5(0.65) 


146 


Outerl 


13 36 10.3 


-29 59 


12.0 


5 


322. 6{ 42.4) 


388.2(16.5) 


31.5(0.28) 


15.52(0.28) 


9.17(0.84) 


21.6(0.39) 


147 


Outerl 


13 36 17.8 


-29 58 


57.1 


8 


1030. 2{ 71.0) 


223.6(27.4) 


26.8(0.47) 


14.98(0.47) 


8.94(1.40) 


8.8(0.65) 


148 


Outerl 


13 36 08.4 


-29 56 


11.9 


6 


561. 8{ 56.7) 


380.8(21.9) 


21.5(0.37) 


12.85(0.37) 


8.66(1.12) 


6.1(0.52) 


149 


Outerl 


13 35 30.3 


-30 02 


12.1 


5 


425. 8{ 49.6) 


342.5(19.2) 


15.2(0.33) 


8.91(0.33) 


5.46(0.98) 


10.2(0.45) 


150 


OuterL 


13 37 16.0 


-30 02 


53.2 


7 


1340.1(118.9) 


2028.5(33.4) 


31.0(0.68) 


26.34(0.66) 


21.16(1.82) 


16.5(1.04) 


151 


OuterL 


13 37 18.1 


-30 05 


13.0 


11 


8140.0(186.4) 


6986.7(52.3) 


176.2(1.06) 


79.08(1.03) 


74.38(2.87) 


113.5(1.63) 


152 


OuterL 


13 37 10.9 


-30 11 


20.1 


7 


610.5(118.6) 


853.9(33.4) 


25.3(0.68) 


18.88(0.66) 


19.36(1.82) 


5.7(1.04) 


Note. 


— (a) Rad 


us of the photometric aperture, 


in units 


of arc second, (b) Fluxes are in 


units of 10~^^ 


ergs cm 


A-'; the 


1 tj error bars 



are indicated in parenthesis, following the flux value. For FUV and NU V band, the flux h as been corrected for the foreground Milky Way extinction, 
E(B-V)^0.066 [Apuv = 8.376E(B - V) and Ajvf yv - 8.741£;(B - V) ^Wvder et al.ll2005H l after convolving the GALEX spectral response curve with 
the Galaxy's extinction curve, using the relation of lCardelli, Clayton fc Mathij il989D with R^— 3.1); the flux in the SpitzcrlRAC bands is assumed to be 
extinction free. 
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4.2. Mass and Age distributions 

For the SED fitting procedure, we adopt the x^-minimization technique to compare the 
observation and the theoretical SEDs of the sources: 



x\t,E{B-V),m,Z) 



{Lobs — A * Ljnodel)'^ 



S ^ 

JV obs 



where N is the number of the filters available for each source (due to the potential dust 
contamination of 5.8 and 8.0 fim bands, the other four bands, FUV, NUV, 3.6 and 4.5 
fim are used, then N=4). Lobs and L^odei are the luminosity from observation and model. 
"A" is the ratio between the actual mass m in each aperture and the default value 10^ 
provided by the Starburst99 models, a obs accounts for the photometric calibration 

and measurement uncertainties. Since we only have four data points available for each 
source, we fix both the metallicity and the mean dust extinction values. Fixing the lat- 
ter, in particular, overcomes the age-extinction degeneracy that plagues stellar continuum 
band fits. We use as a guidance both the HI m ap and the results on HII regions in the 



M83 outer regions from I Gil de Paz et al.l (120071 ). The HI map gives mean column densi 



ties 3.5x10^° cm ^ and 5x10^° cm ^ for Outerl and O uterL, respectively, corresponding 



E(B— V)=0.07 and =0.1 (IBohlin. Savage fc Drakdll978l ). As an extreme case, we also per- 



form fits using E ( B— ^ 1=0.3, which is at the high-end of the extinction values observed by 



Gil de Paz et al.l (120071 ) in these outer regions. High extinction va lues are likely to be ap - 



plicable to the youngest SF regions, of which the HII regions of I Gil de Paz et al.l (120071 ) 
are part, while the lower extinction values are likely to be applicable to older regions. 
Afuv/E{B -V) = 8.376 and Anuv/E{B - V) = 8.741 for the MW type dust have been 
used for calculating the extinction. L^ q and L4.5 are assumed to be dust free. The age and 
mass of each source corresponding to the minimized are listed in table [3l for both the low 
and high extinction values. The best fit values are listed together with the uncertainty range 
corresponding to the 90% confidence level for the best fits; this is derived from requiring that, 
in the age and mass parameter space, Sx'^= 4.61 (90% confidence level for 2 parameters). 
Uncertainties are relatively small for our best fit parameters, owing to the leverage provided 
by the long-wavelength baseline of the data. 

Fig. m shows two examples of the fitting results for each extinction value. The two 
upper panels report the case of a source with a good match between data and models in 
the four primary bands (where a 'good match' is defined as agreement between data and 
models within the 3 a error bar of each data point), while the bottom two panels show the 
case of a source with unusual high excess 4.5, 5.8 and 8.0 fim fiuxes (one of the last 16 
sources in Table [2]). The 68%, 90%, and 99% confidence levels on each best fit (age, mass) 
values are shown in Fig. [5] for two representative regions in our sample. These are both 
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intermediate-luminosity (at 3.6 /iin) sources (see Table |2]), one in Outerl and the other in 
OuterL. The covariance between the two parameters is such that higher masses correspond 
to higher ages. This is readily understood by recalling that the mass is mainly determined 
by the value of the 3.6 fim emission, while the age is constrained by the relative intensity 
of the UV and infrared data; increasing the mass by moving the normalization of the SED 
at higher 3.6 /xm values will then result in a 'flatter' (i.e., older age) best fitting SED to the 
UV data. 




Fig. 4. — Two examples of our SED fitting results. The top two panels show an acceptable fit 
between data and the models, for the two values of the extinction investigated in this work: 
a low value (here E(B— V=0.07 for Outerl) and a high value (E(B— V=0.3). An 'acceptable 
fit' is defined as a match between the models and the data points to within the 3 a error bar 
of each datapoint in the four fitting bands (FUV, NUV, 3.6 fim and 4.5 /im). In this specific 
case, we also observe, as in many other cases, a good match also between the data and the 
models at 5.8 /im, while the 8 /im emission shows a clear excess due to dust emission. The 
bottom two panels show an example of a source with unusual high excess flux at 4.5, 5.8 
and 8.0 /i (one of the last 16 sources in Table [2]). The two sets of lines are the Starburst99 
SEDs without (higher UV values) and with (lower UV values) dust extinction. 
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Table 3. Derived Ages and Masses for the Sources 
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Fig. 5. — Examples of uncertainty diagrams in mass and age for two of our regions, one in 
Outerl (ID69) and the other in OuterL (ID118). The plots are for the low extinction case. 
The lines marking increasing areas show the 68%, 90%, and 99% confidence levels around 
the best fit values. 
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Table 3 — Continued 



Age (Myr) log Mass (Mq) 
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The distributions of the best fit masses and ages for the 152 sources (also separating the 
136 without infrared excess from those with infrared excess, see above) are shown in Fig. [6] 
and Fig. [TJ for both the low and high extinction values. As expected from the age-extinction 
degeneracy, the fits with the higher extinction value produce significantly lower age values, 
as these are mostly determined by the ratio of UV to IR fiux. However, the mean mass, 
which is determined mainly by the 3.6 fim fiux value, is relatively well constrained (to better 
than a factor of 2). For the low extinction values, the mass distribution has a peak around 
10^'^ Mq, with a median value around 10^'^ ^Q- Fo^' high extinction (E(B-V)=0.3), the 
median value decreases to 10^'^ Mq, i.e. by abou t 60% . Masses as large as 10® Mq are 
also found in our sample. In Fig. 3 of iThilker et al.l (120051 ) . the masses of several UV bright 
regions are larger than 10^ Mq, but their derived masses tend to be on average slightly 
smaller than the ones we derive here, possibly owing to the better constraint offered by the 
IRAC measurements. 

Star formation has been an ongoing process in the outskirts of M83 for at least the past 
1 Gyr (Fig. [7]). There is a noticeable increase in the number of sources towards younger 
ages, and the median age for the sample is indeed ~180 Myr; however this increase could 
be due selection effects. Our UV-based source selection method gives preference to young 
stellar population, and we miss most of the old stellar population due to their faint UV 
emission. This effect can be seen clearly from Fig. [HI In this Figure, the distribution 
of age versus mass seems to follow a trend of increasing mean age for increasing mass. 
While the empty region in the bottom-right corner of the plot (regions of old age and low 
total mass) is due to instrumental detection limits, the lack of young and massive stellar 
populations points to possibly a nother type of selection effect, the so-called 'size-of-sample' 
effect (e.g.. iHunter et al.l l2003l . and references therein). For constant star formation, the 
relation between the maximum cluster mass observed and the age of the clusters is related 
by: Mmax oc(age)^/"~^'', where a is the exponent of the initial cluster mass function. By 
fitting a linear relation to the upper envelope of our data in a Log(Mas s)-Log(age) plot, w e 
derive a ~2, a value consistent with results for cluster mass functions (IHunter et al.ll2003l ). 



Within our sample of 136 sources, 13 (~10.5%) have ages of less than 10 Myr. This 
is the age range when ionizing radiation is expected to be produced by massive stars. Our 
observed fraction of 10.5 % ionizing so u rces i s similar to the fraction of UV spots with Ha 
counterparts reported in lThilker et al.l (120051 ). For the higher extinction value (Fig. [71 right 
panel), we obtain that 28.3% sources are younger than 10 Myr, which is mildly inconsistent 
with the dearth of Ha sources in the outer regions of M83. Thus, we infer that E(B-V)~0.1 
is a more realistic extinction value for the majority of the sources in our regions. 



Not surprisingly, the 16 regions with excess 4.5, 5.8, and 8.0 /xm fiuxes have system- 
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Table 3 — Continued 



Age (Myr) log Mass (Mq) 

ID E(B-V)=0.07,0.1'' E(B-V)=0.3'' E(B-V)=0.07,0. 1" E(B-V)=0.3'' 
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20 O + '^ O 


4 44+0-39 


4 24+0-05 
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6 O+l-O 
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3 62+0 03 
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80.0l^S;« 


6 0+1-0 
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136 


140 0+110 
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The sources with unusual high excess 4.5, 5.8 and 8.0 /^m flux 
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^SED fitting results for the 152 sources in the two O uteri and Outer L 
fields. The values of the color excess used in the fits are E(B — V)— 0.07 for 
Outerl and E(B — V) — 0. 1 for Outer L. These values are derived from the aver- 
age HI column densities in each region, iVjf / — 3.5 X 10^*^ cm~^ for Outerl and 
iV^j- — 5 X lO'^'^ cm~^ for OuterL, and the extinction— to— HI column density ratio 
of lBohlin, Savage &: Drakd i 19781 ). 

'^SED fitting results, applying a larger value of the extinction, E(B-V) — 0.3, as 
derived by iGil de Paz et all <200l) for a few of the HII regions detected in the 
MS3 outskirts. 

'^The IDs with added 'XUVn' labels give the cross-identification of the sources 
in common with those of lGil de Paz et al ] i2007tl . 
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Fig. 6. — Histogram of the mass distribution of the sample, as obtained from the SED fitting, 
for both the extinction values adopted: E(B— V)~0.09 (indicative of the actual values of 
E(B-V)~0.07 for Outer! and 0.1 for OuterL, left panel) and E(B-V)=0.3 (right panel. 
Black dash line for the total of 152 sources, blue dot line for the 136 final sources and red 
line for the sources with unusually large 4.5, 5.8 and 8.0 /im excess. 



E(B-V)~0.09 



E(B-V)=0.3 




Fig. 7. — As in Figure [6], but for the age distribution of the sources. 



atically larger masses and on average older ages than the rest of the sample, reinforcing 
the possibility that the IR excess is due to contamination from unrecognized foreground or 
background sources that cannot be separated from the in-situ source of UV emission. In- 
deed, including the 2MASS J,H, and photometry in the SED fitting for these 16 sources 
shows that the 4.5, 5.8, and 8.0 /im IRAC fluxes are systematically larger (in some cases by 
more than an order of magnitude) than the best fit through the GALEX, 2MASS and IRAC 
3.6 fj,m data. Despite the large uncertainties in the 2MASS datapoints, this result confirms 
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Age (Myr) 

Fig. 8. — Mass (in logarithm scale and in units of Mq) as a function of age (in linear scale 
and in units of Myr) for all our 152 sources. The thin lines identified with a bandpass name 
show the 1 a detection limit in that bandpass for a region of 5" radius for the six GALEX 
and Spitzer filters. The actual detection limit for our sources is a complex combination 
of selection criteria involving both UV and infrared detections, as described in section [3l 
In particular, the 3.6 fim detection limit is the dominant selection criterion at young ages, 
while the condition that FUV>5 a dominates at older ages. The 1 a detection limit from 
those combined criteria is shown in the Figure as a thick black line (Final-5") for 5" radius 
apertures and a thick red line (Final-10") for 10" radius apertures. For sources younger than 
~50-100 Myr, the dominant selection criterion is the 3.6 /im detection, while for older ages 
limits in the FUV detections dominate the selection. 

that the last 16 regions of Table [2] are contaminated by (likely) background galaxies at large 
distances. 

The inclusion of the 2MASS datapoints in the SED fitting of the 136 'bonafide' sources 
supports the results obtained for the GALEX+IRAC-only fits: single-age populations are 
still an acceptable fit to the data. The fits that include the 2MASS data produce median 
ages only marginally younger than without the 2MASS data (160 Myr versus 180 Myr), and 
median masses that are a factor 1.6 smaller (10^'''' Mq versus 10^'^ ^o)- '^^^ 2MASS data, 
thus, provide a sanity check that support our baseline results. 

To further test whether even our bonafide 136 final sources may be contaminated by 
unrecognized sources unrelated to M83, we have isolated those sources that are located in 
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regions of HI column densities above the median of the regions (N(HI)>2.2xlO^° cm~^). The 
hypothesis behind this selection is that sources located in correspondence of relatively high 
HI column densities are more likely to be physically associated with M83. There are 60 such 
sources (44% of the total), and their mass and age distributions are shown in Figure O Both 
distributions are consistent with those of the whole sample of 136 sources, with a median 
value for the age which is about half that of the whole sample (~90 Myr versus 180 Myr), 
and a factor 2.5 lower median mass (lO"^'^ Mq versus lO'^'^ Mq). 
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Fig. 9. — Mass (left, in units of Mq) and age (right, in units of Myr) distributions for the 60 
sources located in regions of high HI column density in both Outerl and OuterL. The low 
extinction values have been adopted in the SED fits. 



4.3. Comparison with Ha observations 



Gil de Paz et al.l (120071 ) obtained optical spectra in the Southern parts of the M83 



ou tskirts, and four emissi on line sources fall within our OuterL field: XUVOl-04 (Table 2 
of iGil de Paz et al.l 120071 ). The cross-indentifications with our sources are given in Table [31 
We note significant spatial off-sets between XUV3 and XUV4 and their counterparts in our 
sample (86 and 97, see Table [3]). However, since regions 86 and 97 are the brightest and 
closest UV spots around these two Ha-emitting regions, we believe that they come from the 
same objects. 



From Table [3] and Table HI all regions in common with the sample of iGil de Paz et al. 



( 120071 ) have ages around 100 Myr, which would indicate that they should not be line emitting 



- 26 - 



sources, in apparent contradiction with the results of iGil de Paz et al.l (120071 ) . However, if 
we take as examples XUV2 and XUV3, the extinction values measured by those authors 
for these two regions are 0.23 and 0.29, respectively. If we adopt the larger extinction value 
E(B-V)=0.3 in the SED fitting, we ob tain that XUV2 become s 6 Myr old and XUV3 20 Myr 
old. All four sources in common with iGil de Paz et al.l (120071 ) have indeed ages equal to or 
younger than ~20 Myr for the higher extinction value (Table [3]), which are not incompatible 
with those sources being line-emitting ones. For XUV3, HST images (Thilker et al, 2007, in 
preparation) show the presence of an older (red) stellar population within the source. In our 
IRAC observations, XUV4 is polluted by a nearby foreground star, which may be partially 
responsible for an artificial increase of the deri ved age. XUVl has the l argest H/? fiux among 
the four sources in common with the sample o f iGil de Paz et al.l ( 20071) and has a very small 
extinction as derived from the optical spectra. I Gil de Paz et al.l (120071 ) suggests that the H/5 
fiux of XUVl could be produced by a single young massive star 3 Myr old and with a mass of 
40 Mq. Both the age (80 Myr for the small extinction value, see Table Hj) and the total mass 
(which depends mainly from the IRAC 3.6 /im fiux and implies that the total mass of stars 
above 40 Mq is 4300 Mq) that we derive for XUVl are again in apparent contradiction with 
these results. We suggest that, for XUVl multiple stellar populations, spanning a range of 
ages, are present within the region enclosed by our apertures, with the older ones carrying 
most of the mass and responsible for the IRAC fiux and the younger ones responsible for the 
strong Ha fiux. This hypothesis is also in line with the HST observational results for XUV3. 



4.4. Dust Emission 

Observational evidence for extended dust emission around galaxy disks has increased 
in recent times. Extended (beyond the optical disk) dust emission has been observed in 
a variety of galaxies using both ISO and Spitzer, from large edge-on spirals like NGC891 



Table 4. The counterparts of the Ha sources 



XUV region Namc*^ ID^ fHfUO (ergs s ^ cm *^)^ distance (arcsccond) Agc*^ (Myr) log (Mass^ (solar mass)) 

XUVl 85 1.11x10"^^ 1.35 80. 4.9 

XUV2 90 2.39x10"^'' 4.4 90. 4.7 

XUV3 86 2.07x10^1'^ 12.4 120. 5.4 

XUV4 97 1.60x10^1" 8.1 110. 5.2 



Note. — (a) The XUV region na ming convention from Ta ble 2 of IGil de Paz et al.lt2007D . (b) The ID in table. [2] (c) The H/3 flux 
after the extinction correction from I Gil de Paz et al ] i2007l 1 and (d) The age and mass calculated from our SED fitting method, for 
the low extinction value E{B — V) — 0.1, see Table fsl 
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dPopescu fc Tuffslboosh to dwarfs like UGC10445 jHinz et al.ll2006l l Non-negligible extinc- 
tion has been measured up to 2 effective radi i in spiral galaxies using the overlapping pair 



method (e.g., iHolwerda. Keel fc BoltonI 120071 ) . GALEX observations of the edge-on star- 



burst galaxies M82 and NGC 253 also show unusual high UV luminosity in the halo, which 
cannot be explained by sh ock-heated or photoionized gas and could have dust-scattering 



origi n (IHoopes et al. 



2005r). in agreement with evidence for dust emission in these same ar- 



eas (lEngelbracht et al.l 120061 ). Presence of dust in the outer regions of galaxies can provide 
crucial insights into the metal pollution of the intergalactic medium. 

By calculating the excess infrared flux over the SED best fitting models for our 136 
sources, we get that about 14±6%, 24±30% and 67±26% of the total fiux at 4.5, 5.8 and 8.0 
Hm, respectively, results from the contribution of a non-stellar (dust emission) component. 

The 8 fim dust -only emission can be used to derive an approximate value of the SFR 
in the two regions, keeping in mind a nu r nber of caveats about using the PAH emission 
for tracing star formati on (jCalzetti et al.l 120071 ): it is very sensitive to both metallicity 
(lEngelbracht et al.ll2005l ) and the star formation history of the region under consideration. 
In particular, only when the oxygen abundance is around solar value, there is a relatively 
tight relation between the 8 /im emission and the SFR. Since the abundance in the outer 
disk of M83 is about 1/5-1/10 the solar value, we derive a relation between the 8 una flu x 
and the SFR using the low oxygen abundance data points in Fig. 3 of ICalzetti et al.l (120071 ): 



SFR{Mq yr-') = 1.6 x 10~^'[L8 ^^{erg s-')] 



(2) 



From the 22 sources in both Outerl and OuterL that show an excess at 8 yum (and little 
or no excess in the other bands) we then derive a mean SFR density for the two regions 
by summing up the dust emission at 8 /im and dividing it by the total area. By applying 
equation 2, we obtain an average SFR density of 0.8x10^^ Mq yr~^ kpc~'^. 



(Lu et al.ll2003: 


Helou et al. 


2004 


Reean et al. 



20041 ). This is the flrst time such excess is inferred (albeit with large uncertainty) in the 



outer regions o f a gal axy. The origin of this excess, likely due to dust emission, is still un- 



clear. ILu et al.l (120031 ) suggest that the excess could be due to very small grains transiently 
heated by single photons to high temperatures, ~1000 K. . Our result of 14%±6% by flux as 



due to non-stellar emission is intermediate between the v alue f ound by iRegan et al.l (120041 ) 



for NGC 7 331 ( 6%) and the value found by iHelou et all mm for NGC300 (17%) and by 



Lu et al.l (120031 ) for a sample of galaxies observed with ISO. 
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4.5. The Laws of Star Formation 



As already documented by lBoissier et al.l (120071 ). the steUar and gas radial pr ofiles of the 



galaxies with extended UV emission do obey the scaling laws of star formation (IKennicutt 



19981 : IKennicutt et al.l 120071 ). We investigate the proportionality between the localized star 
formation rate density and gas density in the outer regions of M83, by combining the GALEX 
FUV data and the 8 fim dust-only emission data with the HI map. This analysis thus 
provides a verification of the local scaling laws of star formation in the outskirts of this 
galaxy, in contrast with previous results that have analyzed azimuthal averages of both the 
SFR and gas densities. We will assume that the total SFR in the region is given by the 
sum of the SFRs derived from the observed FUV (uncorrected for extinction) and from the 
8 jUm-du st emission. The 8 /xm-emitting dust is heated by non-ionizing UV and optical 
photons (ILi &: Drainell2002l ). and can thus be adopted as a tracer of the UV photons that 
have been absorbed by dust and re-emitted in the IR; therefore, the extinction-corrected 
UV emission will be the sum of the observed FUV and the 8 /im-dust emission. We will use 
this assumption to derive SFRs in the outskirts of M83. We note that the absence of data on 
the molecular gas content is a limitation of our analysis, which should then be interpreted 
with care. 

At the large distances of Outer! and OuterL from the center of the galaxy, the response 
of the VLA primary beam pattern is substantially reduced (compared to the galaxy center), 
the beam size is degraded, and the intrinsic Nhi column densities are comparatively low. 
To control and minimize uncertainties arising from the combination of these characteristics, 
we perform photometry of regions in apertures of 20" diameter, corresponding to a physical 
scale of about 440 pc, slightly larger than the sizes employed so far in our analysis. We 
use newly defined regions that are measured in both the GALEX FUV and in the dust-only 
8 iim images. The larger apertures enable us to include in each measurement multiple stellar 
clusters, thus justifying the derivation of a SFR density (which would be inappropriate for 
a single stellar cluster). We will also assume that star formation has proceeded at a quasi- 
constant level over the past ^100 Myr, an assumption justified by the results of the previous 
sections. This part of the analysis does not aim at completeness in any sense on the selection 
or measurement of UV- and/or 8 /im-emitting regions; it only aims at providing a range 
of values for the SFRs and HI column densities, to investigate whether star formation at 
such large distances may suggest deviatio ns from the Schmidt-Kenn icutt Law or suggest 
violations of the star formation threshold (IMartin &: Kennicuttll200ll ). We identify a total 
of 54 regions (25 in Outerl and 29 in OuterL), and obtain UV and 8 yum-dust photometry 
for each. 



To infer the critical density for this face-on galaxy, we use the rotation velocity of 
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160 km measured at the distance of ~20 kpc (ISofue et al.lll999l ). assume that the velocity 
dispersion of the stars is neghgib le relative t o the r otation velocity at this distance from the 
center, and apply equation 6 of iKennicutt I (Il989l ). Figure [TD] shows the histogram of the 
ratio of the HI density to the gas critical density at the average distance of the regions 
in Outerl and OuterL. As both fields are located at ~15' ~19.5 kpc distance from the 
center of the galaxy, we adopt a single value of the critical density. The scatter is large, 
as expected for measurements performed at low signal-to-noise ratio, but the peak value is 
around S//7/Scrit=l- This suggests that star formation in these regions happens at roughly 
the local critical density value. The tail towards negative values in Figure [TO] may suggest that 
the HI gas is clumped over scales that are smaller than our measurement aperture (440 pc); 
this is not unreasonable, since the typical HII complex has a much smaller characteristics 
scale than ~400 pc. However, we should stress again that the addition of data on the 
molecular gas density would boost these ratios in the positive sense. 
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Fig. 10. — Histogram of the ratio of the HI gas density to the critical density for selected 
UV-bright regions in the Outerl and OuterL fields. 



The outer regions also display star formation roughly in agreement with the Schmidt- 
Kennicutt Law. This is shown in Figure [TTl where the star formation rate density, SFRD(UV+8 /xm), 
versus HI gas density plot is compa red both with the mean trend observed for RiSOO pc re- 
gio ns in M51 (IKennicutt et al.ll2007r) . and with the average trend for whole galaxi es reported 
in (Kennicutt E"998l ). The UV data are converted to SFRs using the formula of ( Kennicutt 



19981 ) modified for the Kroupa IMF. The conversion between the 8 /im-dust data and SFR 
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uses equation 2. The dynamical range spanned by the M83's outskirts data is small, and 
generally in the low gas density part of the diagram, but there is a general agreement with 
the trend observed for the M51 regions. Yet there is a significant number of M83 regions that 
are located above the mean trend marked by the M51 data, implying that these data show 
higher SFRs for the amount of measured gas density. We should recall that the uncertainties 
in one of the two components of the SFR, i.e., the faint dust-emission-only data, are large 
(Figure 9) and the conversion from 8 fim dust emission to SFR is highly uncertain, due 



to the significant d ependence of the 8 /xm emission on metallicity (lEngelbracht et al.l 12005 



Calzetti et al.ll2007l ). In addition, the M83 data lack measurements of molecular gas content. 



which, if added, would move the data towards higher gas densities. Finally, the obser ved 
trend is not dissimilar from what observed in low-metallicity galaxies (IKennicutt Ill998l ). 




Log[2H] 



Fig. 11. — The star formation rate density (in units of Mq yr^^ kpc^^) versus gas density (in 
units of Mq pc^^) for UV-bright regions (squares) in the Ou terl and OuterL region s, and 
for regions in the disk of M51 (small triangles), as derived in IKennicutt et al.l (120071 ). The 
typical 1 a error bar is shown for the M83 data in the upper-left corner of the plot. The 
gas density for the M83 data is derived from HI alone, while it is a combination of atomic 
and molecular gas measurements for the M51 data. The continuous line is the best fit to 



the M 51 data from IKennicutt et al.l (120071 ): the dashed lines is the relation from IKennicutt 
( 119981 ) for whole galaxies. 
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5. Discussion and Summary 

The new Spitzer IRAC observations presented here, in conjunction with GALEX data 
and a new HI map from THINGS, have expanded our understanding the star formation 
processes in the outskirts of M83. 

The IRAC photometry has targeted the virtually-extinction-free, mid-infrared emission 
from the stars (3.6 and 4.5 /xm), and the dust emission (5.8 and 8.0 /im) from stellar clusters in 
the extreme outer regions of M83, providing constrains on their masses, their ages (together 
with the UV data), and their star formation properties (together with the HI data). 

The Spitzer IRAC observations have targeted two fields, about 9.8 kpc in size at a 
distance of ~19.5 kpc from the center of M83. The two fields are located at about 3 times 
the Ha edge in the galaxy, in correspondence of UV-emitting regions that present a dearth of 
Hq! emission. The cross-identification of UV and IRAC sources has yielded a final sample of 
136 'bona-fide' stellar clusters (or multiple stellar clusters) belonging to M83, after purging 
all potential foreground stellar and background galaxy contaminations with a variety of 
methods. 

Comparison between the multi-wavelength photometry and synthetic SEDs from Star- 
burst99 yields ages for the sources between ~ 1 Myr to more than 1 Gyr , which explains the 



dearth of Ha emission in the outer regions of M83 (iThilker et al.ll2005l ). Our sources show 
a median age around 180 Myr, which however could be a selection effect induced by our 
UV-based source selection technique. Overall, star formation has been an on-going process 
in the outer regions of M83 over the past ~1 Gyr. The SED fits also give masses in the range 
10^ to 10^ Mq for the sources, with a median value of 10^'^ comparable to the masses 
of globular clusters. These results are quantitatively confirmed when the (albeit far more 
uncertain than our IRAC data) 2MASS J, H, and K<j data are added as constraints to the 
SED fitting. It should be remarked that our sources (each encompassing a physical region 
of 220 pc or more in size) could include more than one stellar cluster; this is likely tru e for 



one of the sources in common with the sample of HII regions of I Gil de Paz et al.l (120071 ). for 
which those authors obtain, from optical spectroscopy, masses and ages far lower than what 
we obtain with multi-wavelength SED fitting. 

The agreement between expectations and observations for the ionizing gas emission from 
these regions (i.e., little emission expected overall, due to the broad range of ages shown by 
the stellar clusters) appears to argue against significant leakage of ionizing photons. This 
assumes that there is little dust extinction in those regions. For significant values of the 
extinction, the expected ages of the clusters would decrease, thus introducing a discrepancy 
between the observed and predicted amounts of Ha flux, and the possibility of UV photons 
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leakage. However, measureme nts of the Balmer decrement available for a few of the HII 
regions (IGil de Paz et al.l 120071 ) support our choice of low values for the dust extinction in 
these regions. Presence of ionizing photons leakage would have negligible impact on our 
conclusions on the star formation law in the outer regions, because our measurements of the 
SFR density rely on the non-ionizing UV stellar continuum, both as a direct measure and 
as an indirect measure via the 8 /im dust emission. 

Some of the sources show dust emission at 4.5 /im, 5.8 fim, and 8 /im. We use the 8 /im 
dust emission to derive an average SFR density for the two fields of 0.8 x 10~^ Mq yr~^ kpc'"^., 
although this number should be used with caution, gi ven the many issues related to the use 
of the 8 fim PAH emission for tracing star formation (jCalzetti et al.l 120071 ). The dust emis- 
sion at 4.5 yum we observe accounts for 14%±6% of the total flux; this number is in betwee n 
the values reported for other galaxies (ILu et al.ll2003l : iRegan et al.ll2004l : iHelou et al.l 120041 ). 
although this is th e first time this excess is observed in regions so far removed from the 
center of a galaxy. ILu et al.l (120031 ) suggest that the 4.5 /im excess could be emission from 
very small dust grains transiently heated to 1000 K. The UV photons from the our clusters 
could provide the heating source for the dust, although this will need confirmation. 

Star formation, thus, appears to be an unexceptional event in the outskirts of M83. 
Prom the THINGS HI map the location of our sources is in correspondence of local HI 
enhancements (see also lThilker et al.ll2005l ). and those regions are consistent with the local 
g as density to be around the c ritical density value, in agreement with the threshold hypothesis 
of lMartin fc KennicuttI (120011 ) . at least when applied locally. Furthermore, the star formation 
rate densities and gas densities follo ws a scaling relation s imilar to that found for the star 
forming regions in the disk of M51 (IKennicutt et al.l 120071 ). This reinforces the conclusion 
that the outer regions of M83, albeit sparsely populated with stars and sites of low gas 
densities on average, still form stars and stellar clusters following relations already established 
for the higher density regions of disks. 
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ABSTRACT 



Spitzer IRAC observations of two fields in the extended UV disk (XUV-disk) 
of M83 have been recently obtained, ~3 Rhii away from the center of the galaxy 
{Rhii=Q.6 kpc). GALEX UV images have shown the two fields to host in-situ 
recent star-formation. The IRAC images are used in conjunction with GALEX 
data and new HI imaging from The HI Nearby Galaxy Survey (THINGS) to 
constrain stellar masses and ages of the UV clumps in the fields, and to relate 
the local recent star formation to the reservoir of available gas. Multi-wavelength 
photometry in the UV and mid-IR bands of 136 UV clumps (spatial resolution 
>220 pc) identified in the two target fields, together with model fitting of the 
stellar UV-MIR spectral energy distributions (SED), suggest that the clumps 
cover a range of ages between a few Myr and >1 Gyr with a median value 
around < 100 Myr, and have masses in the range 10^-3x10^ M©, with a peak 
~ 10^'^Mq. The range of observed ages, for which only a small fraction of 
the mass in stars appears to have formed in the past ~10 Myr, agrees with the 
dearth of Ha emission observed in these outer fields. At the location of our IRAC 
fields, the HI map shows localized enhancement and clumping of atomic gas. A 
comparison of the observed star formation with the gas reservoir shows that the 
UV clumps follow the Schmidt-Kennicutt scaling law of star formation, and that 
star formation is occurring in regions with gas densities at approximately (within 
a factor of a few) the critical density value derived according to the Toomre Q 
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gravitational stability criterion. The significant 8 /im excess in several of the 
clumps (16% of the total by number accounting for ~67% of the 8 /xm fiux)) 
provides evidence for the existence o f dust in these remote fields, in agreement 
with results for other galaxies (e.g. iPopescu fc Tuffj l2003l ). Furthermore, we 



observe a relatively small excess of emission at 4.5 /im in the clumps (14%±6% 
by fiux), which suggests contribution from hot small grains (~ lOOOK), as already 
observed in other galaxies. From our data, the outer regions of the M83 galaxy 
disk show evidence of a time-extended star formation history over <1 Gyr, and of 
a moderately chemically-evolved interstellar medium, in agree ment with recent 



findin gs on the metallicity of the outer HII regions of M83 (iGil de Paz et al. 
2007h . 



Subject headings: galaxies: evolution 
ISM — galaxies: star clusters 



galaxies: individual (M83) — galaxies: 



1. Introduction 



Star formation in the outer disks of galaxies has multiple implications for our under- 
standing of the formation and evolution of disks, of the laws that govern star formation, 
and of the interaction of massive stars with the rarefied interstellar medium in those fields. 
Presence of star formation ensures that those external fields are undergoing some chemical 
enrichment. Radiative and mechanical feedback from massive stars may be more efficient in 
a low-density environment, which may thus play a key role in the enrichment of the pristine 



halo. The star formation proce ss in the outer disk takes place late compared to the inner 

Furthermore, the low gas density enables tests of the star 



disk (iMuhoz-Mateos et al 



formation threshold (IMartin fc Kennicuttll200ll ) and of the relation betweeri gas density and 
star formation rate density (the Schmidt-Kennicutt Law, iKennicutt I Il998l ) at the low end 
of the range. Overall, outer disks provide insights into low-density conditions for the star 
formation that may have characterized the early disk formation. 

Deep Ha imaging had already revealed outer disk star formation beyon d two optical 
radii (-R25) in a few nearby galaxies (IFerguson et al.lll998l : iLelievre fc Royll2000l ). Broadband 
obser vations had also shown presence of sig nificant numbers of B stars in the outer disk of 
M31 JCuillandre et allboOlh and NGC 6822 Jde Blok et al.ll2003h . Unlike their counterparts 
in the inner disks, the outer disk HII regions are small, faint and isolated. Thus, the presence 
of suc h HII regions does not repre sent a challenge to the notion that Ha 'edges' exist in most 
disks dMartin fc Kennicuttl[200ll ). 
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More recently, GALEX has shown that the UV pr ofile of disks extends, with a smooth 



trend, well beyond the Ha radius or the R25 of galaxies (IThilker et al.ll2005l : iGil de Paz et al. 



20051 . Thilker et al., 2007, in press, Boissier et al., 2007, in press). XUV-disks are apparently 
common in the present epoch - though not ubiquitous -, as t hey are found at least at a 
low level in about one-quarter of the spiral galaxy population (IZaritsky fc ChristleinI 120071 . 
Thilker et al. 2007, in press). The UV clumps in the outer disk are generally associated with 
extended HI struct ures, and show eyidenc e for metal enrichment, with metallicities in the 



range Z^/S-Zq/IO flGil de 



within extended HI disks fiThilker et al 



'az et al 



2007). Indeed, extended UV disks are generally hosted 



20051 : iGil de Paz et al.ll2005l ). Thus, the rarified 



outer disk is not quiet and star formation process is occurring stealthily in the extended HI 
disk. The presence of smooth UV surface brightness profiles does not counter the paucity of 
Ha emission found in the outer disk UV emission, since the UV probes a larger range of ages 
than Ha and at low-star- formation levels the number of ionizing stars may be very small 
and stochastically absent (IBoissier et al.l 120071 ). Thus, the outer disk star formation merits 
careful investigation, to understand its nature and the conditions under which it occurs. 

The southern grand-design galaxy M83 (RA(2000): 13h37m00.9s, Decf2000): -29d5 1m57s) 



is a metal rich (> Z q) Sc galaxy, viewed almost face on ( i^ 25", ISofue et al.lll999l ). Its 
distance, 4.5 Mpc (IThim et al.l l2003l : iKarachentsevI 120051 ). makes the spatial resolution 
of Spitzer (2", ~44 pc) and GALEX (5", ~109 pc), sufficient to isolate large star for- 
mation complexes. The mean star formation rate per unit area of M83 is ro ughly 0.04 



Mf7\ yr ^ kpc ^ and the total star formation rate (SFR) is about 5 Mq yr ^ (IKennicutt 



19981 ). Therefore, M83 is relatively active among the local normal star forming galaxies. 
The Ha emission of M83 has a wel l-defined 'edge' at the galactocentric distance of around 
6.6 kpc (IMartin fc Kennicuttll200ll ). M83 also has a very extended HI disk which is about 
3 times larger than the optical one (R25). Th e possibihty of a clos e tidal encounter with its 
companion NGC 5253, around 1-2 Gyr ago (IRogstad et al.lll974l : Ivan den Bergh.l Il980l ) is 
still controversial, but if confirmed, it could have triggered the starbursts in both galaxies. 
Additio nal sources of interaction f or this galaxy include dwarf companion galaxies such as 
KK208 jKarachentseva et al.lE998[ ). 



GALEX images of M83 have revealed the presence of more than 100 UV-bright sources 
grouped in highly str uctured complexes beyond the radius where the Ha surface brightness 
decreases abruptly (IThilker et al.ll2005l ). However, only a few of these UV bright regions 
(10%-20%) have Ha counterparts. This apparent discrepancy has been sugges ted to be due 



to th e wide range of the stellar populations' age within the outer disk of M83 (IThilker et al. 



20051 ) and/or to the stochastic nature of the stellar IMF at the very low gas densities in 
these remote regions (IBoissier et al.l 120071 ). 
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New mid-infrared images, from Spitzer IRAC , of two fields in the extreme outskirts of 
M83 are presented here, and are combined with t he GALEX iraages and a 2 1 cm (HI) image 
from The HI Nearby Galaxy Survey (THINGS, IWalter et all J2005I . boosh ) to investigate 
the properties of the UV-bright extended disk; we are exploiting the leverage offered by the 
multi-wavelength observations to constrain ages and masses, and to infer the star formation 
characteristics of the UV knots in the outer regions of M83. The targeted IRAC fields 
show relatively high UV surface brightness and are associated with large scale filamentary 
HI structures. The THINGS HI image offers about 3 times better angular resolution than 
pre-existing HI images, and enables us to explore in detail the scaling laws of star formation 
at a local level. 



In section 2, we provide a brief description of our data. We describe in section 3 the 
methods used to select the sample of regions for the detailed analysis. The analysis of the 
photometry is presented in section 4 and our results are summarized in section 5. 



Observations and Data Reduction 



Spitzer images at 3.6, 4.5, 5.8 and 8.0 /im of two UV bright f ields in the outsk irts of 
M83 were obtained during 2005, July using the IRAC instrument (IFazio et al.ll2004l ). The 
two fields, called Outerl and OuterL, are selected for the combination of both UV-brightness 
and Ha deficiency (see Fig. [1]) and detailed information of these fields is listed in Table [1] 
Spitzer IRAC observed each field twice with a time interval of several days to enable removal 
of asteroids. The two fields were mapped with 4x7 and 4x6 grids, for Outerl and OuterL, 
respectively, with sub-frame steps in order to reach maximum depth in fields of size of 
7.5' X 7.5' and 7.5' x 5'. The total observation time per pixel in the center of each map is 480 
second, reaching a depth of 0.007 MJy/sr at 3.6 /im (1 a detection error). The standard 
Spitzer Infrare d Nearby Galaxies Survey (SINGS) IRAC pipeline was employed to create the 
final mosaics (IKennicutt et al.ll2003l : iRegan et al.ll2004l ). 



M83 was observed by GALEX in 2003 in both the far-UV (FUV, A ~15 29 A) and the 
near-UV (NUV, A ;^231 2 A), as part of the Nearby Galaxy Survery (NGS; iBianchi et al. 



20031 : iMartin et al.ll2005l ). The total exposure time is 1352 s in each filter and the Point 
Spread Function has a FWHM=4.6^^ The detailed description of the GALEX observations 
is provided in iThilker et al.l (120051 ). The GALEX images were aligned, registered, and 
cropped to match the Field-of- Views (FOVs) of the IRAC images of Outerl and OuterL, 
respectively, (Fig. [1]) and re-sampled to the same pixel size (0.75" /pixel for Outerl and 
1" /pixel for OuterL). The GALEX and Spitzer images of Outerl and OuterL are depicted 
in Fig. [2] and Fig. [31 respectively . 



- 5 - 



The 21 cm HI image of M83, from THINGS flWalter et all l2005l . l2008h . was reduced 
and calibrated according to standard pro cedures developed for the THINGS project (e.g. 
Walter et al.l l2005l : iKennicutt et al.l 120071 ). with the addition of a careful blanking in the 
data cube to improve sensitivity in the outer regions of the map. Furthermore, the map has 
been corrected for the primary beam attenuation, as the two IRAC fields lay at the edges of 
the half-power point of VLA primary beam (~30' diameter). The final beamsize is 15"xll". 
The final HI image has also been aligned, registered, and cropped to match both the Outerl 
and Outer! FOVs. 



Images from the 2MASS survey (ISkrutskie et al.ll2006l ) at J, H, and were also used 
in the present analysis to provide a sanity check on our results. The 2MASS images are at 
least 2 orders of magnitude less deep than our IRAC images, and the resulting uncertainty 
on the photometry of our sources (which are very faint in the 2MASS fields) prevents the 
effective use of the 2MASS data in our fitting routines (see below). However, the 2MASS 
data can still provide consistency checks, and they will be used as such. The 2MASS images 
corresponding to the locations of our fields where retrieved and aligned to both Outerl and 
OuterL. 



3. Source Identification and Selection 

Outerl and OuterL contain a few hundred UV sources, which can be, in addition to 
star forming regions in M83, background galaxies and foreground stars. Identification of the 
nature of the sources and of the "bona fide" IRAC counterparts to their UV emission (rather 
than chance superpositions) has required extensive investigation. We developed procedures 
to both exclude distant galaxies and foreground stars (often showing in both the UV and mid- 
IR, M83 lies at an low ecliptic latitude of ~ 18°), and flag UV sources that are contaminated 



Table 1. Spitzer IRAC observations 



Target 


Field Center Position 


Galactocent. 


IRAC Map's 


exposure 


Field 


(J2000) 


Distance 


core size 


tinie(s)'^ 


M83 Outerl 


13:35:56.02 -29:57:23.0 


15.06' (19.71 kpc=2.96 Rhii) 


7.5'x7.5' 


3000 


M83 OuterL 


13:36:56.90 -30:06:41.7 


14.77' (19.33 kpc=2.91 Rhii) 


7.5'x5.0' 


2573 



''For the central map locations, the effective exposure is 480 second 
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Fig. 1. — [adapted from iThilker et al.l (120051)]: The GALEX image of M83 in shown in a 
two-color combination, blue for the FU V channel and orange for the NUV channel. The red 
contours outline the extented HI disk (ITilanus &: Allenlll993l ). The yellow contour indicate 
^neutraigas = 10 Mq pc~^. The two green squares at the right and the bottom of picture indi- 
cate the positions and approximate size of our two Spitzer IRAC fields; the name convention 
for the two fields is also indicated on the figure. North is up, East is left. 
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Fig. 2.— GALEX (FUV, NUV) and Spitzer IRAC (3.6 and 8.0 /xm) images of the Outerl 
field. The black circles show the 97 sources retained in the sample after removal of contami- 
nating foreground and background sources. The size of the circles indicates the apertures in 
which photometry is performed. 
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Fig. 3. — As in Figure [21 now for the OuterL field. A total of 55 sources are retained for 
this field. 

in the IRAC bands by chance superposition of physically unrelated mid-IR sources. We 
applied several steps, as discussed below, to achieve these goals. 
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We began constructing our sample by selecting FUV- and NUV-emitting regions from 
the GALEX images in order to include all potential star forming regions. The FUV is sen- 
sitive to the young phases of SF (though not effectively instantanenous like Ha) , while 
the NUV traces the slightly evolved stell ar population (age ~ 100 Myr). Sextractor (ver- 



sion2.4.4, iBertin. E.. fc Arnouts. S.lll996l ) was employed to identify all 2a sources, yielding 
226 and 133 in Outerl and OuterL, respectively (after manually combining several adjacent 
and overlapping apertures in the crowded regions). Since the instrument background and 
sky background of GALEX and Spitzer are much different, we did not use the Sextractor 
photometry package to measure the flux of these apertures in the IRAC bands. Instead, 
the IRAF routine "phot" was used to extract the photometry of the UV-detected sources in 
all the six GALEX-I- IRAC bands after subtracting the local background from each of them. 
The radius of each region was selected >5" (~the GALEX PSF, corresponding to a physical 
scale ~220 pc in diameter). The same position and aperture diameter was used in all six 
(FUV, NUV, 3.6, 4.5, 5.8 and 8.0 fim) bands and aperture corrections are applied to the re- 
sulting photometry. All our sources appear point-like at the resolution of both GALEX and 
Spitzer/IRAC; therefore we apply to each source/bandpass combination the point-source 
aperture correction appropriate for the chosen region size. For the GALEX measurements, 
we derived the aperture corrections from photometry with increasing aperture radii of point 
sources in the M83 field (outside the galaxy); as a sani t y che ck, we compared our results with 



the aperture corrections published by lMorrissey et al.l (120071 ). and found no major discrepan- 



cies with our values. For the IRAC measurements, we referred to the point-source aperture 



corrections published in the IRAC Data HandbooMj. Some sources are dim in several bands, 
and after subtracting the local background, the flux is negative. For these la upper 

limit has been set. 

Uncertainties on the measurements for each source/bandpass combination are a quadratic 
combination of three contributions: variance of the local background, photometric calibration 
uncertainties, and variations from potential mis-registration of the multiwavelength images. 
We find that in all cases the variance of the local background is the predominant source of 
uncertainty for our measurements, so we only consider this contribution to the uncertainties 
quoted in Table [21 The local background is derived from fits of the pixel distribution of 
regions of a few thousands pixels (or about l/50th of the area of each image) immediately 
surrounding each aperture, for which we derive the mode and the variancqj. Each region 



^IRAC Data Handbook, v3.0, Ch. 5; |http://ssc. spitzer. caltech.edu/irac/dh/ 

^The routine used for the fits, the fRAF task MSKY written by M. Dickinson (1993, private commu- 
nication), aUows the user to interactively define the interval in the pixel distribution where the mode and 
variance are calculated. This ensures that robust results are obtained even in the absence of source masking. 



- 10 - 



was also visually inspected to ensure that no systematic trends in background areas were 
present. In addition, tests were run on sample apertures to verify that the mode of each local 
background is not dependent on the size of the area selected for the background measure- 
ments: we varied region sizes by up a factor of 2 (both increasing and decreasing) without 
encountering for the vast majority of cases variations on excess of the measured variances. 

Contaminating foreground and background sources are eliminated through several meth- 
ods. For some nearby stars, we can identify them by their shapes in the 3.6 /im images, i.e., 
by the presence of diffraction spikes. We also exploit HST ACS images of the center parts 
of the two fields, obtained by one of us (D.T.), to further identify and remove contaminated 
sources from the sample; the ACS images, thanks to their superior angular resolution, en- 
able a quick recognition of foreground stars and background galaxies. For the most part, we 
drop from our candidates list any source with strong contamination. In those cases where 
the contaminating source is physically much smaller than our measurement aperture (e.g., 
a distant background galaxy) and close to its edge, we shift and resize the measurement 
aperture to avoid it. In a few instances, a foreground star may lay close to the center of our 
selected aperture: in these cases, we assume that the FUV and NUV fluxes receive negligible 
contamination, and we only modify our IRAC photometric measurements, by subtracting 
the contaminating contribution measured in an inset small radius aperture. In all cases, we 
apply the point-source aperture correction appropriate for the aperture size selected. 

Three especially bright MW stars lie within the FOVs of Outerl (to the right) and 
OuterL (both right and left), see Fig. [2] and Fig. [31 They have strong wings due to the PSF 
of the Spitzer IRAC instrument. We reject all apertures that fall within and close to these 
wings. After all identified contaminating sources are removed from or corrected for in our 
initial list, there are 166 and 114 regions remaining in the Outerl and OuterL, respectively. 

Finally, in order to remove dim foreground low mass stars more completely, we adopt 
two strong selection criteria: log{j^^) > —0.5 and fpuv > 5o", where fpuv and Jnuv, in 
unit of ergs cm~^ A~^, are the flux of FUV and NUV separately and a is the measurement 
uncertainty. Although the criteria will cause the loss of some clusters with age larger than 
~ 1 Gyr (from the prediction of Starburst99, for an instantaneous burst of star formation 
with 1/5 solar abundance), it is very effective at eliminating low mass stars which emit most 
of their energy in the optical and near-Infrared band. This 'loss' of the oldest clusters will 
have minimal impact on this analysis, which centers on more recent star formation events 
(younger than about 1 Gyr). 

The final size of our sample is 152 sources: 97 in Outerl and 55 in OuterL. The final 
source list with their fluxes and 1 a uncertainties is listed in table [2l Although we try our 
best to exclude the possible foreground stars and background galaxies, there are still about 
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16 sources showing unusual high excess 4.5, 5.8 and 8.0 fim flux which cannot be easily 
explained by stellar population synthesis models (see section These sources tend to 

be among the bright ones in our sample, and we cannot attribute the observed excess to 
photometric or background subtraction uncertainties. We separate these unusual sources 
from the other normal sources in table [2l We consider the remaining 136 sources in our final 
list (although we will carry along the other 16 sources from some of the comparisons, for 
completeness.) 



Table 2. Source List 







RA(J2000) 


Dec(J2000) 




FUV'' 


NUV'' 


3.6^ 


4.5^ 


5.8^ 


8.0'' 


ID 


IRAC field 


degree 


degree 


Radius* 


Fx (lo-) 




Fx 


Fx (!<') 


Fx (l") 


Fx 



1 


Outerl 


13 


36 


22 


1 


-29 


52 


22 


1 


11 


4U4cs.4( yy 


2) 


2700 


2 


38 


5) 


45 


6 





65) 


17 


64 





65) 


7 


27 


1 


96) 


5 
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55 





94) 


2 


54 


2 


60) 





3 


1 


48) 


97 


OutcrL 


13 


36 


56 


7 


-30 


06 


48 


2 


10 


7707.0 
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4. Analysis 



4.1. Stellar populations' modelling 

The age and mass of each source is constrained by fitting the UV+MIR bands with 
models of the spectral energy distribution (SEDs) of single-age stellar populations. This 
implicitly assumes that the emission in the IRAC bands is due to pho tospheric eniission 



from stars. This is a reasonable assumption for the 3.6 and 4.5 /im bands (jPahre et al.l 12004 



Calzetti et al.ll2005l ). The 5.8 and 8 fim band potentially contain a significant (or dominant) 
contribution from dust emission, due to Policyclic Aromatic Hydrocarbons single-photon 
heated by the stellar populations; to evaluate the impact of dust emission on our stellar 
population fits, we create a dust-only 8 /xm image by rescahng and subtracting the 3.6 /im 
image from the 8 fim one. The resulting dust-only image shows a number of sources with 
detected dust emission, in most cases a weak detection, in agreement wi th the expectation 



that low-metallicity regio ns have underluminous 8 fim dust emission ( lEngelbracht et al. 



20051 : ICalzetti et al.l 120071 ). However, to prevent out fitting results from being even mildly 



contaminated by dust emission contributions, we use only the two IRAC bands at the shortest 
wavelengths in our fits, together with the two GALEX bands (next section). 

The theoretical SEDs are frorn the Starburst99 stellar population synthesis models 
( jLeitherer et al.lll999l : IVdzquez et al.ll2005l ). We select instantaneous burst populations, since 
we are likely to deal with single or small groups of stellar clusters (i.e., same or similar age 
populations) due to the small physical sizes of the regions sampled (lOOpc to 300pc). A 
metalhcity value of Z=0.004 is adopted for the models, to match as closely as possible the 
oxy gen abundances of roug hly 1/10-1/5 Zq measured in a few of the HIl regions in these ar- 
eas (iGil de Paz et al.ll2007l ). The Kroupa Initial Mass Function ($( m)(xm~°, with a=1 .3 for 
O.IMq < M < 0.5Mq, and a=2.3 for 0.5Mo < M < lOOM©, i JKroupa et allboOlh . the 
default IMF of Starburst99, is employed. Since we include in the fits the mid- Infrared (MIR) 
fluxes of each stellar complex, we use the new Padova stellar models with full asymptotic 
giant branch evolution. In order to balance the accuracy and the speed of the simulations for 
different age ranges, we select different time steps for the models: 0.05 Myr for 0.1-1 Myr, 
0.1 Myr for 1-10 Myr, 0.5 Myr for 10-50 Myr, 1 Myr for 50-200 Myr, 2 Myr for 200-500 Myr, 
10 Myr for 500-1000 Myr and 50 Myr for 1-2 Gyr. We convolve the spectral response curve 
of the GALEX and Spitzer filter bandpasses with the simulated stellar population spectra 
to obtain the synthetic photometry in the six bands. Finally, the luminosities and the colors 
from observations are compared with the theoretical model to constrain ages and masses. 
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Table 2 — Continued 







RA{J2000) 


Dec(J2000) 




FUV'' 


NUV'' 


3.6'' 


4.5'' 


5.S^ 


8.0*= 


ID 


IRAC field 


degree 


degree 


Radius^ 


Fx (l"^) 




Fx {1<^) 


Fx (1<^) 


Fx (l"^) 


Fx (1<t) 


131 


OuterL 


13 37 14.8 


-30 09 


55.7 


7 


773.0(118.7) 


546.5(33.4) 


5.1(0.68) 


0.67(0.66) 


0.92(1.82) 


0.3(1.04) 


132 


OuterL 


13 36 56.5 


-30 07 31.0 


6 


1022.4(101.8) 


630.9(28.6) 


0.4(0.58) 


0.21(0.57) 


0.67(1.56) 


0.3(0.89) 


133 


OuterL 


13 36 58.0 


-30 07 


52.9 


5 


823. 0( 84.8) 


509.5(23.9) 


1.2(0.48) 


0.56(0.47) 


0.41(1.30) 


0.5(0.74) 


134 


OuterL 


13 36 57.9 


-30 08 


20.8 


5 


635. 2( 84.8) 


486.8(23.9) 


2.8(0.48) 


1.22(0.47) 


1.44(1.30) 


0.4(0.74) 


135 


OuterL 


13 36 57.1 


-30 07 


59.0 


6 


536.5(101.8) 


873.6(28.6) 


65.8(0.58) 


27.69(0.57) 


9.00(1.56) 


2.5(0.89) 


136 


OuterL 


13 36 57.1 


-30 09 


19.8 


5 


665. 3( 84.8) 


352.5(23.9) 


3.5(0.48) 


1.81(0.47) 


0.89(1.30) 


0.8(0.74) 










Sources with unusual high excess 


4.5, 5.8 and 8 


flux 








137 


Outerl 


13 36 03.4 


-29 55 


25.7 


8 


1420. 2( 70.8) 


742.6(27.4) 


57.5(0.47) 


29.98(0.47) 


15.08(1.39) 


39.5(0.65) 


138 


Outerl 


13 35 55.1 


-29 56 


57.5 


7 


3702. 0{ 63.9) 


1937.4(24.7) 


13.7(0.42) 


7.54(0.42) 


8.18(1.26) 


5.7(0.58) 


139 


Outerl 


13 36 13.4 


-29 59 


42.1 


8 


1474. 3{ 71.0) 


2028.5(27.4) 


61.2(0.47) 


61.33(0.47) 


51.61(1.40) 


35.7(0.65) 


140 


Outerl 


13 35 49.9 


-30 00 


40.7 


8 


5022. 6{ 71.0) 


5299.1(27.4) 


56.0(0.47) 


55.45(0.47) 


50.67(1.40) 


40.0(0.65) 


141 


Outerl 


13 35 28.4 


-29 58 


39.7 


8 


1907.3{ 71.0) 


1437.9(27.4) 


52.0(0.47) 


21.95(0.47) 


13.67(1.40) 


24.9(0.65) 


142 


Outerl 


13 35 39.8 


-29 58 


48.3 


8 


1405. 0{ 71.0) 


1066.6(27.4) 


75.6(0.47) 


43.50(0.47) 


19.21(1.40) 


65.1(0.65) 


143 


Outerl 


13 36 15.6 


-29 58 


00.0 


8 


493. 6{ 71.0) 


864.1(27.4) 


93.3(0.47) 


60.74(0.47) 


41.36(1.40) 


68.2(0.65) 


144 


Outerl 


13 35 48.9 


-29 56 


40.9 


8 
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145 
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8 
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61.5(0.47) 


34.68(0.47) 


18.38(1-40) 


20.5(0.65) 


146 
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12.0 


5 
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388.2(16.5) 


31.5(0.28) 


15.52(0.28) 


9.17(0.84) 


21.6(0.39) 


147 
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57.1 


8 
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223.6(27.4) 


26.8(0.47) 


14.98(0.47) 


8.94(1.40) 


8.8(0.65) 


148 


Outerl 
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11.9 


6 


561. 8{ 56.7) 


380.8(21.9) 


21.5(0.37) 


12.85(0.37) 


8.66(1.12) 


6.1(0.52) 


149 


Outerl 


13 35 30.3 


-30 02 


12.1 


5 


425. 8{ 49.6) 


342.5(19.2) 


15.2(0.33) 


8.91(0.33) 


5.46(0.98) 


10.2(0.45) 


150 


OuterL 


13 37 16.0 


-30 02 


53.2 


7 


1340.1(118.9) 


2028.5(33.4) 


31.0(0.68) 


26.34(0.66) 


21.16(1.82) 


16.5(1.04) 


151 


OuterL 


13 37 18.1 


-30 05 


13.0 


11 


8140.0(186.4) 


6986.7(52.3) 


176.2(1.06) 


79.08(1.03) 


74.38(2.87) 


113.5(1.63) 


152 


OuterL 


13 37 10.9 


-30 11 


20.1 


7 


610.5(118.6) 


853.9(33.4) 


25.3(0.68) 


18.88(0.66) 


19.36(1.82) 


5.7(1.04) 


Note. 


— (a) Rad 


us of the photometric aperture, 


in units 


of arc second, (b) Fluxes are in 


units of 10~^^ 


ergs cm 


A-'; the 


1 tj error bars 



are indicated in parenthesis, following the flux value. For FUV and NU V band, the flux h as been corrected for the foreground Milky Way extinction, 
E(B-V)^0.066 [Apuv = 8.376E(B - V) and Ajvf yv - 8.741£;(B - V) ^Wvder et al.ll2005H l after convolving the GALEX spectral response curve with 
the Galaxy's extinction curve, using the relation of lCardelli, Clayton fc Mathij il989D with R^— 3.1); the flux in the SpitzcrlRAC bands is assumed to be 
extinction free. 
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4.2. Mass and Age distributions 

For the SED fitting procedure, we adopt the x^-minimization technique to compare the 
observation and the theoretical SEDs of the sources: 



x\t,E{B-V),m,Z) 



{Lobs — A * Ljnodel)'^ 



S ^ 

JV obs 



where N is the number of the filters available for each source (due to the potential dust 
contamination of 5.8 and 8.0 fim bands, the other four bands, FUV, NUV, 3.6 and 4.5 
fim are used, then N=4). Lobs and L^odei are the luminosity from observation and model. 
"A" is the ratio between the actual mass m in each aperture and the default value 10^ 
provided by the Starburst99 models, a obs accounts for the photometric calibration 

and measurement uncertainties. Since we only have four data points available for each 
source, we fix both the metallicity and the mean dust extinction values. Fixing the lat- 
ter, in particular, overcomes the age-extinction degeneracy that plagues stellar continuum 
band fits. We use as a guidance both the HI m ap and the results on HII regions in the 



M83 outer regions from I Gil de Paz et al.l (120071 ). The HI map gives mean column densi 



ties 3.5x10^° cm ^ and 5x10^° cm ^ for Outerl and O uterL, respectively, corresponding 



E(B— V)=0.07 and =0.1 (IBohlin. Savage fc Drakdll978l ). As an extreme case, we also per- 



form fits using E ( B— ^ 1=0.3, which is at the high-end of the extinction values observed by 



Gil de Paz et al.l (120071 ) in these outer regions. High extinction va lues are likely to be ap - 



plicable to the youngest SF regions, of which the HII regions of I Gil de Paz et al.l (120071 ) 
are part, while the lower extinction values are likely to be applicable to older regions. 
Afuv/E{B -V) = 8.376 and Anuv/E{B - V) = 8.741 for the MW type dust have been 
used for calculating the extinction. L^ q and L4.5 are assumed to be dust free. The age and 
mass of each source corresponding to the minimized are listed in table [3l for both the low 
and high extinction values. The best fit values are listed together with the uncertainty range 
corresponding to the 90% confidence level for the best fits; this is derived from requiring that, 
in the age and mass parameter space, Sx'^= 4.61 (90% confidence level for 2 parameters). 
Uncertainties are relatively small for our best fit parameters, owing to the leverage provided 
by the long-wavelength baseline of the data. 

Fig. m shows two examples of the fitting results for each extinction value. The two 
upper panels report the case of a source with a good match between data and models in 
the four primary bands (where a 'good match' is defined as agreement between data and 
models within the 3 a error bar of each data point), while the bottom two panels show the 
case of a source with unusual high excess 4.5, 5.8 and 8.0 fim fiuxes (one of the last 16 
sources in Table [2]). The 68%, 90%, and 99% confidence levels on each best fit (age, mass) 
values are shown in Fig. [5] for two representative regions in our sample. These are both 
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intermediate-luminosity (at 3.6 /iin) sources (see Table |2]), one in Outerl and the other in 
OuterL. The covariance between the two parameters is such that higher masses correspond 
to higher ages. This is readily understood by recalling that the mass is mainly determined 
by the value of the 3.6 fim emission, while the age is constrained by the relative intensity 
of the UV and infrared data; increasing the mass by moving the normalization of the SED 
at higher 3.6 /xm values will then result in a 'flatter' (i.e., older age) best fitting SED to the 
UV data. 




Fig. 4. — Two examples of our SED fitting results. The top two panels show an acceptable fit 
between data and the models, for the two values of the extinction investigated in this work: 
a low value (here E(B— V=0.07 for Outerl) and a high value (E(B— V=0.3). An 'acceptable 
fit' is defined as a match between the models and the data points to within the 3 a error bar 
of each datapoint in the four fitting bands (FUV, NUV, 3.6 fim and 4.5 /im). In this specific 
case, we also observe, as in many other cases, a good match also between the data and the 
models at 5.8 /im, while the 8 /im emission shows a clear excess due to dust emission. The 
bottom two panels show an example of a source with unusual high excess flux at 4.5, 5.8 
and 8.0 /i (one of the last 16 sources in Table [2]). The two sets of lines are the Starburst99 
SEDs without (higher UV values) and with (lower UV values) dust extinction. 
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Table 3. Derived Ages and Masses for the Sources 



Age (Myr) 
E(B-V)=0. 07,0.1" B(B-V)=0.3'' 



log Mass (Mq) 
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23 
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28 
29 
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31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
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Fig. 5. — Examples of uncertainty diagrams in mass and age for two of our regions, one in 
Outerl (ID69) and the other in OuterL (ID118). The plots are for the low extinction case. 
The lines marking increasing areas show the 68%, 90%, and 99% confidence levels around 
the best fit values. 
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Table 3 — Continued 



Age (Myr) log Mass (Mq) 
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The distributions of the best fit masses and ages for the 152 sources (also separating the 
136 without infrared excess from those with infrared excess, see above) are shown in Fig. [6] 
and Fig. [TJ for both the low and high extinction values. As expected from the age-extinction 
degeneracy, the fits with the higher extinction value produce significantly lower age values, 
as these are mostly determined by the ratio of UV to IR fiux. However, the mean mass, 
which is determined mainly by the 3.6 fim fiux value, is relatively well constrained (to better 
than a factor of 2). For the low extinction values, the mass distribution has a peak around 
10^'^ Mq, with a median value around 10^'^ ^Q- Fo^' high extinction (E(B-V)=0.3), the 
median value decreases to 10^'^ Mq, i.e. by abou t 60% . Masses as large as 10® Mq are 
also found in our sample. In Fig. 3 of iThilker et al.l (120051 ) . the masses of several UV bright 
regions are larger than 10^ Mq, but their derived masses tend to be on average slightly 
smaller than the ones we derive here, possibly owing to the better constraint offered by the 
IRAC measurements. 

Star formation has been an ongoing process in the outskirts of M83 for at least the past 
1 Gyr (Fig. [7]). There is a noticeable increase in the number of sources towards younger 
ages, and the median age for the sample is indeed ~180 Myr; however this increase could 
be due selection effects. Our UV-based source selection method gives preference to young 
stellar population, and we miss most of the old stellar population due to their faint UV 
emission. This effect can be seen clearly from Fig. [HI In this Figure, the distribution 
of age versus mass seems to follow a trend of increasing mean age for increasing mass. 
While the empty region in the bottom-right corner of the plot (regions of old age and low 
total mass) is due to instrumental detection limits, the lack of young and massive stellar 
populations points to possibly a nother type of selection effect, the so-called 'size-of-sample' 
effect (e.g.. iHunter et al.l l2003l . and references therein). For constant star formation, the 
relation between the maximum cluster mass observed and the age of the clusters is related 
by: Mmax oc(age)^/"~^'', where a is the exponent of the initial cluster mass function. By 
fitting a linear relation to the upper envelope of our data in a Log(Mas s)-Log(age) plot, w e 
derive a ~2, a value consistent with results for cluster mass functions (IHunter et al.ll2003l ). 



Within our sample of 136 sources, 13 (~10.5%) have ages of less than 10 Myr. This 
is the age range when ionizing radiation is expected to be produced by massive stars. Our 
observed fraction of 10.5 % ionizing so u rces i s similar to the fraction of UV spots with Ha 
counterparts reported in lThilker et al.l (120051 ). For the higher extinction value (Fig. [71 right 
panel), we obtain that 28.3% sources are younger than 10 Myr, which is mildly inconsistent 
with the dearth of Ha sources in the outer regions of M83. Thus, we infer that E(B-V)~0.1 
is a more realistic extinction value for the majority of the sources in our regions. 



Not surprisingly, the 16 regions with excess 4.5, 5.8, and 8.0 /xm fiuxes have system- 
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Table 3 — Continued 



Age (Myr) log Mass (Mq) 

ID E(B-V)=0.07,0.1'' E(B-V)=0.3'' E(B-V)=0.07,0. 1" E(B-V)=0.3'' 
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^SED fitting results for the 152 sources in the two O uteri and Outer L 
fields. The values of the color excess used in the fits are E(B — V)— 0.07 for 
Outerl and E(B — V) — 0. 1 for Outer L. These values are derived from the aver- 
age HI column densities in each region, iVjf / — 3.5 X 10^*^ cm~^ for Outerl and 
iV^j- — 5 X lO'^'^ cm~^ for OuterL, and the extinction— to— HI column density ratio 
of lBohlin, Savage &: Drakd i 19781 ). 

'^SED fitting results, applying a larger value of the extinction, E(B-V) — 0.3, as 
derived by iGil de Paz et all <200l) for a few of the HII regions detected in the 
MS3 outskirts. 

'^The IDs with added 'XUVn' labels give the cross-identification of the sources 
in common with those of lGil de Paz et al ] i2007tl . 
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Fig. 6. — Histogram of the mass distribution of the sample, as obtained from the SED fitting, 
for both the extinction values adopted: E(B— V)~0.09 (indicative of the actual values of 
E(B-V)~0.07 for Outer! and 0.1 for OuterL, left panel) and E(B-V)=0.3 (right panel. 
Black dash line for the total of 152 sources, blue dot line for the 136 final sources and red 
line for the sources with unusually large 4.5, 5.8 and 8.0 /im excess. 



E(B-V)~0.09 



E(B-V)=0.3 




Fig. 7. — As in Figure [6], but for the age distribution of the sources. 



atically larger masses and on average older ages than the rest of the sample, reinforcing 
the possibility that the IR excess is due to contamination from unrecognized foreground or 
background sources that cannot be separated from the in-situ source of UV emission. In- 
deed, including the 2MASS J,H, and photometry in the SED fitting for these 16 sources 
shows that the 4.5, 5.8, and 8.0 /im IRAC fluxes are systematically larger (in some cases by 
more than an order of magnitude) than the best fit through the GALEX, 2MASS and IRAC 
3.6 fj,m data. Despite the large uncertainties in the 2MASS datapoints, this result confirms 
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Fig. 8. — Mass (in logarithm scale and in units of Mq) as a function of age (in linear scale 
and in units of Myr) for all our 152 sources. The thin lines identified with a bandpass name 
show the 1 a detection limit in that bandpass for a region of 5" radius for the six GALEX 
and Spitzer filters. The actual detection limit for our sources is a complex combination 
of selection criteria involving both UV and infrared detections, as described in section [3l 
In particular, the 3.6 fim detection limit is the dominant selection criterion at young ages, 
while the condition that FUV>5 a dominates at older ages. The 1 a detection limit from 
those combined criteria is shown in the Figure as a thick black line (Final-5") for 5" radius 
apertures and a thick red line (Final-10") for 10" radius apertures. For sources younger than 
~50-100 Myr, the dominant selection criterion is the 3.6 /im detection, while for older ages 
limits in the FUV detections dominate the selection. 

that the last 16 regions of Table [2] are contaminated by (likely) background galaxies at large 
distances. 

The inclusion of the 2MASS datapoints in the SED fitting of the 136 'bonafide' sources 
supports the results obtained for the GALEX+IRAC-only fits: single-age populations are 
still an acceptable fit to the data. The fits that include the 2MASS data produce median 
ages only marginally younger than without the 2MASS data (160 Myr versus 180 Myr), and 
median masses that are a factor 1.6 smaller (10^'''' Mq versus 10^'^ ^o)- '^^^ 2MASS data, 
thus, provide a sanity check that support our baseline results. 

To further test whether even our bonafide 136 final sources may be contaminated by 
unrecognized sources unrelated to M83, we have isolated those sources that are located in 
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regions of HI column densities above the median of the regions (N(HI)>2.2xlO^° cm~^). The 
hypothesis behind this selection is that sources located in correspondence of relatively high 
HI column densities are more likely to be physically associated with M83. There are 60 such 
sources (44% of the total), and their mass and age distributions are shown in Figure O Both 
distributions are consistent with those of the whole sample of 136 sources, with a median 
value for the age which is about half that of the whole sample (~90 Myr versus 180 Myr), 
and a factor 2.5 lower median mass (lO"^'^ Mq versus lO'^'^ Mq). 
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Fig. 9. — Mass (left, in units of Mq) and age (right, in units of Myr) distributions for the 60 
sources located in regions of high HI column density in both Outerl and OuterL. The low 
extinction values have been adopted in the SED fits. 



4.3. Comparison with Ha observations 



Gil de Paz et al.l (120071 ) obtained optical spectra in the Southern parts of the M83 



ou tskirts, and four emissi on line sources fall within our OuterL field: XUVOl-04 (Table 2 
of iGil de Paz et al.l 120071 ). The cross-indentifications with our sources are given in Table [31 
We note significant spatial off-sets between XUV3 and XUV4 and their counterparts in our 
sample (86 and 97, see Table [3]). However, since regions 86 and 97 are the brightest and 
closest UV spots around these two Ha-emitting regions, we believe that they come from the 
same objects. 



From Table [3] and Table HI all regions in common with the sample of iGil de Paz et al. 



( 120071 ) have ages around 100 Myr, which would indicate that they should not be line emitting 
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sources, in apparent contradiction with the results of iGil de Paz et al.l (120071 ) . However, if 
we take as examples XUV2 and XUV3, the extinction values measured by those authors 
for these two regions are 0.23 and 0.29, respectively. If we adopt the larger extinction value 
E(B-V)=0.3 in the SED fitting, we ob tain that XUV2 become s 6 Myr old and XUV3 20 Myr 
old. All four sources in common with iGil de Paz et al.l (120071 ) have indeed ages equal to or 
younger than ~20 Myr for the higher extinction value (Table [3]), which are not incompatible 
with those sources being line-emitting ones. For XUV3, HST images (Thilker et al, 2007, in 
preparation) show the presence of an older (red) stellar population within the source. In our 
IRAC observations, XUV4 is polluted by a nearby foreground star, which may be partially 
responsible for an artificial increase of the deri ved age. XUVl has the l argest H/? fiux among 
the four sources in common with the sample o f iGil de Paz et al.l ( 20071) and has a very small 
extinction as derived from the optical spectra. I Gil de Paz et al.l (120071 ) suggests that the H/5 
fiux of XUVl could be produced by a single young massive star 3 Myr old and with a mass of 
40 Mq. Both the age (80 Myr for the small extinction value, see Table Hj) and the total mass 
(which depends mainly from the IRAC 3.6 /im fiux and implies that the total mass of stars 
above 40 Mq is 4300 Mq) that we derive for XUVl are again in apparent contradiction with 
these results. We suggest that, for XUVl multiple stellar populations, spanning a range of 
ages, are present within the region enclosed by our apertures, with the older ones carrying 
most of the mass and responsible for the IRAC fiux and the younger ones responsible for the 
strong Ha fiux. This hypothesis is also in line with the HST observational results for XUV3. 



4.4. Dust Emission 

Observational evidence for extended dust emission around galaxy disks has increased 
in recent times. Extended (beyond the optical disk) dust emission has been observed in 
a variety of galaxies using both ISO and Spitzer, from large edge-on spirals like NGC891 



Table 4. The counterparts of the Ha sources 



XUV region Namc*^ ID^ fHfUO (ergs s ^ cm *^)^ distance (arcsccond) Agc*^ (Myr) log (Mass^ (solar mass)) 

XUVl 85 1.11x10"^^ 1.35 80. 4.9 

XUV2 90 2.39x10"^'' 4.4 90. 4.7 

XUV3 86 2.07x10^1'^ 12.4 120. 5.4 

XUV4 97 1.60x10^1" 8.1 110. 5.2 



Note. — (a) The XUV region na ming convention from Ta ble 2 of IGil de Paz et al.lt2007D . (b) The ID in table. [2] (c) The H/3 flux 
after the extinction correction from I Gil de Paz et al ] i2007l 1 and (d) The age and mass calculated from our SED fitting method, for 
the low extinction value E{B — V) — 0.1, see Table fsl 
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dPopescu fc Tuffslboosh to dwarfs like UGC10445 jHinz et al.ll2006l l Non-negligible extinc- 
tion has been measured up to 2 effective radi i in spiral galaxies using the overlapping pair 



method (e.g., iHolwerda. Keel fc BoltonI 120071 ) . GALEX observations of the edge-on star- 



burst galaxies M82 and NGC 253 also show unusual high UV luminosity in the halo, which 
cannot be explained by sh ock-heated or photoionized gas and could have dust-scattering 



origi n (IHoopes et al. 



2005r). in agreement with evidence for dust emission in these same ar- 



eas (lEngelbracht et al.l 120061 ). Presence of dust in the outer regions of galaxies can provide 
crucial insights into the metal pollution of the intergalactic medium. 

By calculating the excess infrared flux over the SED best fitting models for our 136 
sources, we get that about 14±6%, 24±30% and 67±26% of the total fiux at 4.5, 5.8 and 8.0 
Hm, respectively, results from the contribution of a non-stellar (dust emission) component. 

The 8 fim dust -only emission can be used to derive an approximate value of the SFR 
in the two regions, keeping in mind a nu r nber of caveats about using the PAH emission 
for tracing star formati on (jCalzetti et al.l 120071 ): it is very sensitive to both metallicity 
(lEngelbracht et al.ll2005l ) and the star formation history of the region under consideration. 
In particular, only when the oxygen abundance is around solar value, there is a relatively 
tight relation between the 8 /im emission and the SFR. Since the abundance in the outer 
disk of M83 is about 1/5-1/10 the solar value, we derive a relation between the 8 una flu x 
and the SFR using the low oxygen abundance data points in Fig. 3 of ICalzetti et al.l (120071 ): 



SFR{Mq yr-') = 1.6 x 10~^'[L8 ^^{erg s-')] 



(2) 



From the 22 sources in both Outerl and OuterL that show an excess at 8 yum (and little 
or no excess in the other bands) we then derive a mean SFR density for the two regions 
by summing up the dust emission at 8 /im and dividing it by the total area. By applying 
equation 2, we obtain an average SFR density of 0.8x10^^ Mq yr~^ kpc~'^. 



(Lu et al.ll2003: 


Helou et al. 


2004 


Reean et al. 



20041 ). This is the flrst time such excess is inferred (albeit with large uncertainty) in the 



outer regions o f a gal axy. The origin of this excess, likely due to dust emission, is still un- 



clear. ILu et al.l (120031 ) suggest that the excess could be due to very small grains transiently 
heated by single photons to high temperatures, ~1000 K. . Our result of 14%±6% by flux as 



due to non-stellar emission is intermediate between the v alue f ound by iRegan et al.l (120041 ) 



for NGC 7 331 ( 6%) and the value found by iHelou et all mm for NGC300 (17%) and by 



Lu et al.l (120031 ) for a sample of galaxies observed with ISO. 
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4.5. The Laws of Star Formation 



As already documented by lBoissier et al.l (120071 ). the steUar and gas radial pr ofiles of the 



galaxies with extended UV emission do obey the scaling laws of star formation (IKennicutt 



19981 : IKennicutt et al.l 120071 ). We investigate the proportionality between the localized star 
formation rate density and gas density in the outer regions of M83, by combining the GALEX 
FUV data and the 8 fim dust-only emission data with the HI map. This analysis thus 
provides a verification of the local scaling laws of star formation in the outskirts of this 
galaxy, in contrast with previous results that have analyzed azimuthal averages of both the 
SFR and gas densities. We will assume that the total SFR in the region is given by the 
sum of the SFRs derived from the observed FUV (uncorrected for extinction) and from the 
8 jUm-du st emission. The 8 /xm-emitting dust is heated by non-ionizing UV and optical 
photons (ILi &: Drainell2002l ). and can thus be adopted as a tracer of the UV photons that 
have been absorbed by dust and re-emitted in the IR; therefore, the extinction-corrected 
UV emission will be the sum of the observed FUV and the 8 /im-dust emission. We will use 
this assumption to derive SFRs in the outskirts of M83. We note that the absence of data on 
the molecular gas content is a limitation of our analysis, which should then be interpreted 
with care. 

At the large distances of Outer! and OuterL from the center of the galaxy, the response 
of the VLA primary beam pattern is substantially reduced (compared to the galaxy center), 
the beam size is degraded, and the intrinsic Nhi column densities are comparatively low. 
To control and minimize uncertainties arising from the combination of these characteristics, 
we perform photometry of regions in apertures of 20" diameter, corresponding to a physical 
scale of about 440 pc, slightly larger than the sizes employed so far in our analysis. We 
use newly defined regions that are measured in both the GALEX FUV and in the dust-only 
8 iim images. The larger apertures enable us to include in each measurement multiple stellar 
clusters, thus justifying the derivation of a SFR density (which would be inappropriate for 
a single stellar cluster). We will also assume that star formation has proceeded at a quasi- 
constant level over the past ^100 Myr, an assumption justified by the results of the previous 
sections. This part of the analysis does not aim at completeness in any sense on the selection 
or measurement of UV- and/or 8 /im-emitting regions; it only aims at providing a range 
of values for the SFRs and HI column densities, to investigate whether star formation at 
such large distances may suggest deviatio ns from the Schmidt-Kenn icutt Law or suggest 
violations of the star formation threshold (IMartin &: Kennicuttll200ll ). We identify a total 
of 54 regions (25 in Outerl and 29 in OuterL), and obtain UV and 8 yum-dust photometry 
for each. 



To infer the critical density for this face-on galaxy, we use the rotation velocity of 



- 29 - 



160 km measured at the distance of ~20 kpc (ISofue et al.lll999l ). assume that the velocity 
dispersion of the stars is neghgib le relative t o the r otation velocity at this distance from the 
center, and apply equation 6 of iKennicutt I (Il989l ). Figure [TD] shows the histogram of the 
ratio of the HI density to the gas critical density at the average distance of the regions 
in Outerl and OuterL. As both fields are located at ~15' ~19.5 kpc distance from the 
center of the galaxy, we adopt a single value of the critical density. The scatter is large, 
as expected for measurements performed at low signal-to-noise ratio, but the peak value is 
around S//7/Scrit=l- This suggests that star formation in these regions happens at roughly 
the local critical density value. The tail towards negative values in Figure [TO] may suggest that 
the HI gas is clumped over scales that are smaller than our measurement aperture (440 pc); 
this is not unreasonable, since the typical HII complex has a much smaller characteristics 
scale than ~400 pc. However, we should stress again that the addition of data on the 
molecular gas density would boost these ratios in the positive sense. 
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Fig. 10. — Histogram of the ratio of the HI gas density to the critical density for selected 
UV-bright regions in the Outerl and OuterL fields. 



The outer regions also display star formation roughly in agreement with the Schmidt- 
Kennicutt Law. This is shown in Figure [TTl where the star formation rate density, SFRD(UV+8 /xm), 
versus HI gas density plot is compa red both with the mean trend observed for RiSOO pc re- 
gio ns in M51 (IKennicutt et al.ll2007r) . and with the average trend for whole galaxi es reported 
in (Kennicutt E"998l ). The UV data are converted to SFRs using the formula of ( Kennicutt 



19981 ) modified for the Kroupa IMF. The conversion between the 8 /im-dust data and SFR 
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uses equation 2. The dynamical range spanned by the M83's outskirts data is small, and 
generally in the low gas density part of the diagram, but there is a general agreement with 
the trend observed for the M51 regions. Yet there is a significant number of M83 regions that 
are located above the mean trend marked by the M51 data, implying that these data show 
higher SFRs for the amount of measured gas density. We should recall that the uncertainties 
in one of the two components of the SFR, i.e., the faint dust-emission-only data, are large 
(Figure 9) and the conversion from 8 fim dust emission to SFR is highly uncertain, due 



to the significant d ependence of the 8 /xm emission on metallicity (lEngelbracht et al.l 12005 



Calzetti et al.ll2007l ). In addition, the M83 data lack measurements of molecular gas content. 



which, if added, would move the data towards higher gas densities. Finally, the obser ved 
trend is not dissimilar from what observed in low-metallicity galaxies (IKennicutt Ill998l ). 




Log[2H] 



Fig. 11. — The star formation rate density (in units of Mq yr^^ kpc^^) versus gas density (in 
units of Mq pc^^) for UV-bright regions (squares) in the Ou terl and OuterL region s, and 
for regions in the disk of M51 (small triangles), as derived in IKennicutt et al.l (120071 ). The 
typical 1 a error bar is shown for the M83 data in the upper-left corner of the plot. The 
gas density for the M83 data is derived from HI alone, while it is a combination of atomic 
and molecular gas measurements for the M51 data. The continuous line is the best fit to 



the M 51 data from IKennicutt et al.l (120071 ): the dashed lines is the relation from IKennicutt 
( 119981 ) for whole galaxies. 
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5. Discussion and Summary 

The new Spitzer IRAC observations presented here, in conjunction with GALEX data 
and a new HI map from THINGS, have expanded our understanding the star formation 
processes in the outskirts of M83. 

The IRAC photometry has targeted the virtually-extinction-free, mid-infrared emission 
from the stars (3.6 and 4.5 /xm), and the dust emission (5.8 and 8.0 /im) from stellar clusters in 
the extreme outer regions of M83, providing constrains on their masses, their ages (together 
with the UV data), and their star formation properties (together with the HI data). 

The Spitzer IRAC observations have targeted two fields, about 9.8 kpc in size at a 
distance of ~19.5 kpc from the center of M83. The two fields are located at about 3 times 
the Ha edge in the galaxy, in correspondence of UV-emitting regions that present a dearth of 
Hq! emission. The cross-identification of UV and IRAC sources has yielded a final sample of 
136 'bona-fide' stellar clusters (or multiple stellar clusters) belonging to M83, after purging 
all potential foreground stellar and background galaxy contaminations with a variety of 
methods. 

Comparison between the multi-wavelength photometry and synthetic SEDs from Star- 
burst99 yields ages for the sources between ~ 1 Myr to more than 1 Gyr , which explains the 



dearth of Ha emission in the outer regions of M83 (iThilker et al.ll2005l ). Our sources show 
a median age around 180 Myr, which however could be a selection effect induced by our 
UV-based source selection technique. Overall, star formation has been an on-going process 
in the outer regions of M83 over the past ~1 Gyr. The SED fits also give masses in the range 
10^ to 10^ Mq for the sources, with a median value of 10^'^ comparable to the masses 
of globular clusters. These results are quantitatively confirmed when the (albeit far more 
uncertain than our IRAC data) 2MASS J, H, and K<j data are added as constraints to the 
SED fitting. It should be remarked that our sources (each encompassing a physical region 
of 220 pc or more in size) could include more than one stellar cluster; this is likely tru e for 



one of the sources in common with the sample of HII regions of I Gil de Paz et al.l (120071 ). for 
which those authors obtain, from optical spectroscopy, masses and ages far lower than what 
we obtain with multi-wavelength SED fitting. 

The agreement between expectations and observations for the ionizing gas emission from 
these regions (i.e., little emission expected overall, due to the broad range of ages shown by 
the stellar clusters) appears to argue against significant leakage of ionizing photons. This 
assumes that there is little dust extinction in those regions. For significant values of the 
extinction, the expected ages of the clusters would decrease, thus introducing a discrepancy 
between the observed and predicted amounts of Ha flux, and the possibility of UV photons 
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leakage. However, measureme nts of the Balmer decrement available for a few of the HII 
regions (IGil de Paz et al.l 120071 ) support our choice of low values for the dust extinction in 
these regions. Presence of ionizing photons leakage would have negligible impact on our 
conclusions on the star formation law in the outer regions, because our measurements of the 
SFR density rely on the non-ionizing UV stellar continuum, both as a direct measure and 
as an indirect measure via the 8 /im dust emission. 

Some of the sources show dust emission at 4.5 /im, 5.8 fim, and 8 /im. We use the 8 /im 
dust emission to derive an average SFR density for the two fields of 0.8 x 10~^ Mq yr~^ kpc'"^., 
although this number should be used with caution, gi ven the many issues related to the use 
of the 8 fim PAH emission for tracing star formation (jCalzetti et al.l 120071 ). The dust emis- 
sion at 4.5 yum we observe accounts for 14%±6% of the total flux; this number is in betwee n 
the values reported for other galaxies (ILu et al.ll2003l : iRegan et al.ll2004l : iHelou et al.l 120041 ). 
although this is th e first time this excess is observed in regions so far removed from the 
center of a galaxy. ILu et al.l (120031 ) suggest that the 4.5 /im excess could be emission from 
very small dust grains transiently heated to 1000 K. The UV photons from the our clusters 
could provide the heating source for the dust, although this will need confirmation. 

Star formation, thus, appears to be an unexceptional event in the outskirts of M83. 
Prom the THINGS HI map the location of our sources is in correspondence of local HI 
enhancements (see also lThilker et al.ll2005l ). and those regions are consistent with the local 
g as density to be around the c ritical density value, in agreement with the threshold hypothesis 
of lMartin fc KennicuttI (120011 ) . at least when applied locally. Furthermore, the star formation 
rate densities and gas densities follo ws a scaling relation s imilar to that found for the star 
forming regions in the disk of M51 (IKennicutt et al.l 120071 ). This reinforces the conclusion 
that the outer regions of M83, albeit sparsely populated with stars and sites of low gas 
densities on average, still form stars and stellar clusters following relations already established 
for the higher density regions of disks. 
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